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EXPLOSIVE  LOADING  OF  METALS 


I.  Preface 

The  explosive  loading  of  metals,  especially  as  related  to  shaped-charge 
liners  and  metals  driven  by  explosives,  is  a  vast  field.  Thousands  of 
documents  have  been  written  covering  many  aspects  of  the  explosive  loading  of 
materials.  However,  due  to  the  intense  interest  in  this  field  by  the  defense 
establishments  of  many  nations,  a  great  deal  of  the  literature  is  either 
classified,  not  available  in  English,  or  hidden  in  obscure  government  or 
industry  documents.  Also,  oertain  Industries  (e.g.,  drilling,  mining,  metal 
forming,  etc.)  engaged  in  explosive  loading,  do  not  readily  distribute 
information  due  to  their  highly  competitive  nature.  This  increases  the  number 
of  proprietary  and  minimal  distribution  documents.  In  short,  a  great  deal  of 
information  is  simply  not  available  to,  or  cannot  be  found  by,  even  the  more 
diligent  researcher. 

I  am  certainly  not  aware  of  all  the  explosive  loading  information 
available.  Furthermore,  I  will  not  even  attempt  to  present  all  the 
information  at  my  disposal.  Instead,  I  will  present  basic  introductory 
material  and  cover  only  the  major  aspects  of  the  explosive  loading  of  metals 
technology  and  limit  my  discussion  to  shaped-charges  and  the  motion  of  metals 
driven  by  explosives.  The  material  and  literature  covered  will  be  that  of 
important  historical  value  and  that  material  found  to  be  most  useful  to  me. 

My  sources  stem  primarily  from  my  association  with  Ballistic  Research 
Laboratory.  As  a  direct  consequence  of  this  association,  I  am  probably  remiss 
in  reporting  the  results  of  other  government  and  industrial  agencies, 
especially  those  which  do  not  deal  directly  with  the  Ballistic  Research 
Laboratory.  Also,  reference  to  foreign  literature  will  be  minimized.  Several 
sources  have  been  omitted  as  being  beyond  the  scope  of  this  text  or  merely 
redundant  or  supplemental  to  material  already  included.  Nevertheless,  a  large 
shaped-charge  formation  and  penetration  bibliography  is  given.  Again,  I  do 
not  claim  that  this  bibliography  is  complete. 

Finally,  a  note  of  thanks  to  Lisa  Ann  Weismiller  for  the  meticulous 
typing  of  a  difficult  manuscript. 


II.  The  Shaped-Charge  Concept 

A  cylinder  of  explosive  with  a  hollow  cavity  in  one  end  and  a  detonator 
at  the  opposite  end  is  known  as  a  hollow  charge.  The  hollow  cavity,  which  may 
assume  almost  any  geometric  shape  such  as  a  hemisphere,  a  cone  or  the  like 
causes  the  gaseous  products  formed  from  the  initiation  of  the  explosive  at  the 
end  of  the  cylinder  opposite  the  hollow  cavity  to  focus  the  energy  of  the 
detonation  products.  The  focusing  of  the  detonation  products  creates  an 
intense  localized  force.  This  concentrated  force  when  directed  against  a 
metal  plate  is  capable  of  creating  a  deeper  cavity  than  a  cylinder  of 
explosive  without  a  hollow  cavity,  even  though  more  explosive  is  available  in 
the  latter  case.  As  discussed  later,  this  phenomenon  is  known  in  the  U.S.  as 
the  Munroe  effect.  Further,  if  the  hollow  cavity  is  lined  with  a  thin  layer 
of  metal,  glass,  ceramic  or  the  like,  the  liner  forms  a  jet  when  the  explosive 
charge  is  detonated.  This  focusing  of  the  liner  material  causes  a  still 
deeper  cavity  to  be  formed  when  the  explosive  is  detonated  with  the  base  of 
the  liner  in  contact  with  a  metal  plate.  The  cavity  formed  becomes  deeper  yet 
when  the  explosive  charge  containing  the  liner  is  removed  some  distance  away 
from  the  plate.  This  distance,  for  which  an  optimum  exists,  is  called  the 
standoff  distance.  Devices  of  this  nature  are  called  lined  cavity  charges  or 
shaped-charges. 

In  the  following  sections,  these  concepts  will  be  pursued  in  detail. 
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III.  History  of  Shaped-Charges 

This  section  will  discuss  shaped-charge  jets  or  the  shaped-charge  effect 
resulting  from  the  interaction  of  an  explosive  with  a  metallic  liner. 

Further,  the  discussion  will  be  restricted  to  conventional  shaped-charge 
rounds,  e.g.,  conical,  hemispherical  and  other  arcuate  metallic  liners,  and 
Self  Forging  Fragments  (SFF)  or  Explosively  Formed  Penetrators  (EFP)  or  P 
charges  will  only  be  briefly  addressed.  Explosive  welding  and  cutting  or 
tapping  charges  will  not  be  considered  in  any  detail. 

The  history  of  shaped-charge  development  is  somewhat  ambiguous  in  that 
the  British,  Germans  and  Americans  all  have  made  significant  claims  to  the 
early  development  of  lined  cavity  charges.  I  understand  that  currently  Dr. 
Manfred  Held  of  MBB,  West  Germany  and  H.  H.  Mohaupt  of  the  U.S.  are  both 
preparing  manuscripts  on  the  history  and/or  development  of  Shaped-Charge 
research.  Dr.  P.  C.  Chou  [1]  and  D.  F.  Kennedy  [2]  have  both  drafted 
manuscripts  on  their  version  of  the  history  and/or  development  of 
shaped-charges.  The  individuals  mentioned  above  are  senior  experts  in  the 
field  and  will  certainly  provide  an  excellent  account  of  current  and/or 
historical  shaped-charge  developments.  Also,  Eather  and  Griffiths  of  the 
li.K.  [3]  have  provided  a  history  of  the  U.K,  contribution  to  the  field  of 
shaped-charges  which  includes  the  contributions  of  Evans,  Ubbelohde,  Taylor, 
Tuck,  Mott,  Hill,  Pack  and  others  which  I  will  discuss  later.  Dr.  Held  may 
include  the  discussions  of  References  1,  2  and  3  in  his  manuscript. 

Recently,  D.  F.  Kennedy  published  a  description  of  the  shaped-charge 
concept  [4]  including  jet  formation,  a  brief  history  of  the  shaped-charge 
concept,  jet  disturbance  mechanisms  and  modern  day  applications  of 
shaped-charge  devices.  I  highly  recommend  Kennedy  [4]  as  required  reading  for 
the  novice  in  the  shaped-charge  field. 

We  will  define  the  term  shaped-charge  to  apply  to  explosive  charges  with 
lined  or  unlined  cavities.  The  cavity  is  formed  in  the  end  of  the  explosive 
charge  opposite  the  point  of  initiation.  The  term  shaped-charge,  however,  has 
a  more  general  meaning,  e.g.,  Cook  [53.  The  shaped-charge  is  also  referred  to 
as  the  hollow  charge  (in  the  U.K.),  the  cumulative  charge  (in  the  USSR),  or 
the  Hohlladung  (in  Germany). 

The  hollow  cavity  effect  was  first  observed  by  Max  von  Foerster  in 
1683  [6,73  and  apparently  rediscovered  by  Charles  Munroe  in  1888 
[8,9.10,11,123.  The  hollow  charge  or  cavity  effect  is  known  in  the  U.S.  and 
U.K.  as  the  Munroe  effect  after  Munroe's  well  documented  experiments  in  1688. 
Munroe  detonated  blocks  of  explosive  in  contact  with  steel  plates.  The 
explosive  charge  had  the  initials  U.S.N.  (United  States  Navy)  inscribed  on  the 
charge  opposite  the  point  of  detonation.  These  Initials  were  reproduced  on 
the  steel  plate.  Munroe  further  observed  that  when  a  cavity  was  formed  in  a 
block  of  explosive,  opposite  the  point  of  initiation,  the  penetration,  or 
depth  of  the  crater  produced  in  the  target,  increased.  In  other  words,  a 
deeper  cavity  could  be  produced  in  a  target  using  a  smaller  mass  of  explosivel 
The  increase  in  penetration  results  from  the  focusing  of  the  explosive  gases 
(detonation  products)  by  the  hollow  cavity.  The  depth  of  the  crater  in  the 
target  can  be  further  increased  by  displacing  the  hollow  charge  some  optimal 
distance  from  the  target,  i.e.,  Increasing  the  standoff-distance  for  a  lined 
cavity  charge. 
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Early  German  reference  to  the  hollow  cavity  effect,  excluding  von 
Foerster,  occurred  in  1911  patents  in  Germany  and  the  U.K,  [13,14].  Also,  in 
1886,  Gustov  Bloem  [15]  of  Dusseldorf  patented  a  shell  for  detonating  caps 
(U.S.  Patent  No.  342,  423)  which  resembles  a  shaped-charge  with  a 
hemispherical  liner. 

Early  British  development  of  the  hollow  cavity  warhead  was  reported  in 
1913  according  to  Reference  [16].  The  earliest  known  Soviet  mention  of  the 
shaped-charge  effect  was  by  Sucharewski  in  1925  [17].  The  first  Italian  paper 
on  the  shaped-charge  effect  was  by  Lodati  in  1932  [18].  Kennedy  [2]  provides 
an  excellent  early  history  of  the  shaped-charge  effect  including  evidence  that 
this  effect  may  have  been  known  prior  to  von  Foerster  using  low  explosives  in 
mining  applications  such  as  mudcapping. 

Studies  involving  the  shaped-charge  effect  continued  from  1911  to  the 
late  1930’s,  mostly  by  German  scientists,  see,  e.g.,  Kennedy  [2].  The 
contributions  of  Watson  [19]  on  percussion  fuzes  and  Wood  [20]  on  SFF,  EFP, 
Misznay-Schardin  devices,  ballistic  discs  or  P-charge  projectiles  appeared 
during  this  period. 

The  lined-cavity  shaped-charge  research  accelerated  tremendously  between 
1935  and  1950  due  primarily  to  World  War  II  and  the  application  of 
shaped-charges  to  the  bazooka,  panzer faust  and  other  devices.  The  discoverers 
of  the  lined  cavity  effect  were  Franz  Rudolf  Thomanek  for  Germany  and  the 
Swiss,  Henry  Hans  Mohaupt  for  the  U.K.  and  the  U.S.  (Unless  we  credit  Bloem 
[15]  with  this  discovery).  At  any  rate,  Thomanek  and  Mohaupt  perfected  this 
concept  and  developed  the  first  effective  lined  cavity  shaped-charge 
penetrators.  Dr.  Thomanek's  early  work  is  given  in  References  [21,22,23,24] 
and  Dr.  Mohaupt' s  work  is  given  in  References  [25,26,27,28].  Kennedy  provides 
additional  references  and  discussion  relating  to  the  work  of  Mohaupt  and 
Thomanek  [2].  Kennedy  also  points  out  that  Munroe  used  a  form  of  a  lined 
cavity  charge  consisting  of  an  open  tin  can  surrounded  by  dynamite  sticks  to 
pierce  a  steel  safe  [12]. 

Mohaupt,  using  lined  cavity  charges,  designed  practical  military  devices 
ranging  from  rifle  grenades,  to  mortars,  to  100  mm  diameter  artillery 
projectiles.  These  devices  were  test  fired  at  the  Swiss  Army  Proving  Ground 
at  Thun,  at  Mohaupt' s  Laboratory  and  at  the  French  Naval  Artillery  Proving 
Ground  at  Gavre.  These  results  were  also  demonstrated  to  the  U.K.  who  then 
began  development  programs  of  their  own.  Following  the  early  results  of  World 
War  II,  the  French  Government  authorized  the  release  of  Mohaupt' s  information 
to  the  U.S.  and  in  late  1940  tests  were  conducted  at  Aberdeen  Proving  Ground, 
Maryland  using  several  aspects  of  lined  cavity  shaped-charges  [28].  The  U.S. 
accepted  the  program,  classified  it,  and  thus  excluded  Dr.  Mohaupt  from  the 
effort  but  produced  the  2.36  in  HEAT  machine  gun  grenade  and  the  75  mm  and  105 
mm  HEAT  artillery  projectiles  in  1941.  Later  the  machine  gun  grenade  was 
modified  to  include  a  rocket  motor  and  a  shoulder  launcher  and  became  the 
bazooka.  The  bazooka  was  first  used  by  the  U.K.  in  North  Africa  in  1941. 

Other  HEAT  rounds  were  fired  from  tank  mounted  howitzers  [2,283. 

The  military  research,  introduced  by  Mohaupt,  was  supported  by  private 
firms  such  as  DuPont  and  the  Carnegie  Institute  of  Technology.  The  Carnegie 
group  (C.I.T.)  employed  some  outstanding  researchers  which  contributed  much  of 
the  current  shaped-charge  knowledge.  The  leaders  at  Carnegie  were 
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Heine-Geldern,  N.  Rostoker,  Emerson  Pugh  and  his  student  Robert  Eichelberger. 
Eichelberger  is  the  former  Director  of  the  Ballistic  Research  Laboratory. 

Shaped-charge  theory  continued  to  develop  during  the  1950's  boosted  by 
the  Korean  War  [5,29,30,31,32,33,34,35,363.  During  this  time  period 
tremendous  progress  was  made  toward  the  understanding  of  the  phenomena 
associated  with  shaped-charge  jets.  Photographic  techniques  were  employed  to 
observe  the  jet  process  and  analytical  models  were  developed.  Efforts  were 
made  to  improve  existing  shaped-charge  liners;  to  use  detonation  wave  shapers; 
to  provide  spin  compensation  via  fluted  liners;  to  provide  shaped-charge 
follow-through  mechanisms;  and  to  enhance  the  overall  system  performance. 

Currently,  shaped-charge  research  continues  in  order  to  devise  successful 
countermeasures  to  the  advanced  armors  currently  fielded  and/or  contemplated, 
see  e.g.,  Kennedy  [43.  Studies  which  originated  in  the  1950's  still  continue; 
notably,  torpedo  applications  of  shaped-charge  rounds,  anti-aircraft  rounds, 
fragmentation  rounds,  multi-staged  or  tandem  liners,  long  standoff  rounds, 
non-conical  liners  and  non-copper  liners.  Also  metallurgical  and  chemical 
aspects  of  the  liner  material  was  well  as  methods  of  liner  fabrication  remain 
important. 

Starting  in  the  1950-60's  significant  shaped-charge  developments  were 
made  possible  by  the  advent  of  excellent  experimental  techniques  such  as  high 
speed  photography  and  flash  radiography.  Other  improvements  resulted  from  the 
transition  from  TNT  to  more  energetic  explosives,  i.e.,  from  TNT  to  Comp  B  to 
Octol  and  then  to  pressed  explosives,  notably  LX-14.  Also,  alternate  modes  of 
initiation  (other  than  point-initiation)  and  waveshaping  techniques  have 
provided  warhead  design  improvements.  Other  advanoes  stemmed  from  the 
development  of  large  computer  codes  to  simulate  the  collapse,  formation  and 
growth  of  the  jet  from  a  shaped-charge  liner.  These  codes  provide  for  the 
most  part,  excellent  descriptions  of  the  formation  of  the  jet.  The  codes, 
however,  are  somewhat  hindered  by  the  lack  of  accurate  equations  of  state  of 
the  liner  material  and  the  explosive  and  by  inadequate  constitutive 
relationships  (i.e.,  the  relationship  betweens  stress,  strain,  strain  rate, 
and  the  velocity  and  velocity  gradients) . 

The  following  sections  will  address  the  reoent  developments  in  more 
detail  and  will  attempt  to  explain  the  shaped— charge  principles  and  related 
areas.  The  information  and  references  given  thus  far  provide  a  brief 
background  of  shaped-charge  development.  Shaped— charge  utilization  was  made 
possible  by  the  discovery  of  blasting  caps  (detonators)  by  Alfred  Nobel 
[28,373  in  1867.  The  explosive  reaction  initiated  by  these  blasting  caps 
oould  propagate  through  a  column  of  explosive  without  the  use  of  confinement. 
This  was  termed  detonation  or  brisant  explosion  [283.  The  hollow  charge 
effect,  described  earlier,  was  exploited  around  1910  by  von  Neumann  £283  and 
was  developed  extensively  during  World  War  II. 

The  Watson  patents  [193  greatly  enhanced  the  weaponlzatlon  and  further 
applications  of  the  shaped-charge  principle.  The  Watson  percussion  fuze 
consisted  of  a  parabolic  booster  charge  with  a  metal  lined  hemispherical  Cor 
arched  shield")  cavity  at  the  "output"  end  of  the  fuze  body  to  intensify  the 
effect  of  the  booster  charge  even  over  an  airgap.  This  fuze  is,  in  effect,  a 
detonator  using  the  shaped-charge  principle. 
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Further  developments  followed  by  Wood  [20]  who  described  what  is  known 
today  as  a  self  forging  fragment,  explosively  formed  penetrator,  a  P-charge  or 
a  Misznay-Schardin  device. 

The  shaped-charge  principle  was  clarified  and  understood  as  a  result  of 
the  pioneering  flash  x-ray  photographs  taken  by  L,  B.  Seely  and  J.  C.  Clark  of 
Ballistic  Research  Laboratory  and  J.  L.  Tuck  of  England,  in  1943,  Based  on  an 
analysis  of  these  photographs,  analytical  models  of  the  collapse  of  a  lined 
conical  shaped-charge  were  developed  by  Birkhoff  [38],  Pugh  [39],  Evans  [40] 
and  others  as  indicated  in  References  [1  and  5],  Important  contributions  were 
carried  out  by  the  Carnegie  group  mentioned  previously  and  continued  by  the 
Ballistic  Research  Laboratory  (L.  Zernow,  R.  Eichelberger  and  Associates). 
Other  laboratories  making  important  contributions  during  this  time  period  were 
the  Naval  Ordnance  Laboratory,  Maryland  (Solem  and  August),  the  Naval  Ordnance 
Test  Station,  California  (Throner,  Weinland,  Kennedy,  Pearson  and  Rinehart), 
Picatinny  Arsenal,  New  Jersey  (Dunkle),  the  Standford  Research  Institute, 
California  (Poulter)  and  others.  The  liner  collapse  and  jet  formation  will  be 
discussed  in  detail  in  a  later  section. 

A  bibliography  and  account  of  the  weaponization  of  the  shaped-charge  and 
similar  principles  is  given  by  J.  Backofen  [41-47],  Backofen’s  bibliography 
is  extensive,  especially  regarding  foreign  sources.  Earlier,  WW  II  results 
and  bibliographical  information  are  given  in  References  [48,49,50], 
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IV.  Introduction  to  Shaped-Charges 

The  directional  penetration  effect  observed  when  a  hollow  charge  is 
detonated  in  contact  with  a  steel  plate  is  graphically  depicted  in  Figure  1 . 
The  crater  depth  is  about  one-half  of  the  diameter  of  the  hollow  of  the 
conical  cavity.  The  cavity  is  produced  by  high  pressure,  high  velocity  gas 
erosion  (the  Munroe  effect).  When  the  hollow  cavity  is  lined  with  a  thin 
hollow  metallic  or  glass  cone,  the  lined  charge  results  in  a  much  deeper 
crater.  Furthermore,  when  the  lined  cavity  charge  is  displaced  from  the 
target  block  some  distance  (known  as  the  standoff)  the  penetration  increases 
even  more.  Figure  1  depicts  these  three  cases. 

The  increase  in  penetration  resulting  from  the  lined  shaped-charge  is  due 
to  the  plastic  deformation  of  the  liner  material  under  the  intense  pressures 
and  temperatures  generated  by  the  detonation  of  the  high  explosive.  The 
mechanism  of  jet  formation  for  metallic  conical  liners  with  a  semi-angle 

(one-half  of  the  conical  apex  angle)  less  than  60°  is  as  described  below. 

Figure  2  shows  a  typical  shaped-charge  warhead  like  that  described  above. 
Note  that  the  explosive  charge  is  not  cylindrical,  but  tapered.  This  removal 
of  some  of  the  explosive  weight  is  termed  "boattailing"  and  does  not  affect 
the  jet  collapse  mechanism.  We  need  only  that  the  detonation  wave  front  be 
plane  and  perpendicular  to  the  longitudinal  axis-of-symmetry. 

When  the  detonator  is  fired,  the  detonation  wave  propagates  through  the 
explosive  with  the  detonation  velocity  of  the  particular  explosive  used.  When 
the  detonation  front  reaches  the  conical  liner,  the  liner  is  subjected  to  the 
intense  pressures  and  high  temperatures  of  the  front  and  begins  to  collapse. 
The  collapse  is  depicted  in  Figure  3  for  the  conical  lined  shaped-charge  shown 
in  Figure  2.  For  the  position  of  the  detonation  wave  front  shown  in  Figure  3, 
the  upper  (apex)  region  of  the  cone  has  collapsed  and  collided  on  the 
axis-of-symmetry .  This  collision  results  in  plastically  deformed  liner 
material  under  tremendous  pressure  being  extruded  along  the  axis-of-symmetry 
[4, HO, 513.  This  extruded  material  is  called  the  jet.  When  the  pressures 
generated  exceed  the  yield  strength  of  the  liner  material,  the  liner  behaves 
approximately  as  an  inviscid,  Incompressible  fluid.  The  cone  collapses 
progressively  from  apex  to  base  under  point  initiation  of  the  high  explosive. 

A  portion  of  the  liner  flows  Into  a  compact  slug  which  is  the  large ,  massive 
portion  at  the  rear  of  the  Jet. 

A  rough  analogy  can  be  drawn  to  the  effect  produced  when  a  sphere  is 
dropped  into  water.  During  impact  and  downward  motion  in  the  water,  the 
cavity  formed  moves  outwards  around  the  sphere  and  then  reverses  to  collide 
along  the  axis-of-symmetry.  On  collision,  vertical  jets  are  formed  which 
squirt  upwards  and  downwards  [40,52). 

Crude  "rules  of  thumb"  were  established  by  Evans  [40]  in  1950.  Namely, 
about  20$  of  the  inner  surface  of  the  metal  of  the  cones  is  used  to  form  the 
jet.  The  Jet  diameter  is  about  one  twentieth  of  the  diameter  of  the  cone. 

The  jet  tip  velocity  is  of  the  order  of  the  detonation  velocity  of  the 
explosive.  The  jet  velocity  decay  is  a  nearly  liner  function  down  to  about 
one  quarter  of  the  detonation  velocity  at  the  tail  (rear)  of  the  Jet.  The 
velocity  of  the  slug  is  of  the  order  of  one  tenth  of  the  Jet  tip  velocity. 
These  estimates  are  only  crude  order-of-magnitude  estimates  and  are  only  for 
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moderate  apex  angle,  copper,  conical  liners.  As  we  will  see,  there  are 
techniques  available  to  provide  more  accurate  estimates  of  the  jet  and  slug 
parameters.  Also  note  that  the  Jet  formation  is  strongly  dependent  on  the 
liner  geometry,  liner  material,  high  explosive  geometry,  confinement  geometry 
and  material  (if  confinement  is  present),  the  type  of  high  explosive  used  and 
the  mode  of  initiation.  These  factors  will  be  addressed  later.  In  any  case, 
the  goal  is  to  direct  and  concentrate  energy  in  the  axial  direction  to  enhance 
the  damage  resulting  from  the  hollow  charge. 

For  any  liner  design,  a  proper  match  between  the  charge  to  mass  ratio 
(explosive  charge  mass  to  liner  mass  ratio)  is  critical.  If  the  liner  is  too 
thick,  the  energy  losses  resulting  from  internal  friction  and  heating  of  the 
liner  walls  during  the  collapse  and  the  energy  losses  due  to  spallation  of  the 
thick  liner  will  reduce  the  collision  velocity  below  the  value  necessary  for 
jet  extrusion.  Also,  if  the  liner  wall  thickness  is  too  thin,  directed  flow 
is  not  achieved  due  to  the  loss  of  structural  integrity  of  the  liner.  If  the 
wall  is  extremely  thin,  the  liner  material  may  undergo  vaporization  upon 
collision. 

Shaped-charges  with  wide  angle  cones  or  hemispherical  liners  show  a 
radically  different  collapse  pattern.  Hemispherical  liners  invert  (or  turn 
inside  out)  from  the  pole  and  as  the  detonation  wave  progresses  toward  the 
base  of  the  liner,  the  hemispherical  liner  approximates  a  conical  liner  and 
the  inverted  collapse  pattern  reverts  to  that  of  a  cone.  In  general,  the 
hemispherical  and  large  angle  conical  liners  usually  result  in  larger 
diameter,  but  lower  velocity  gradient  Jets.  No  massive  slug,  per  se,  is 
formed. 

Of  course,  Jet  effects  or  controlled  fragmentation  are  not  limited  to 
conical  or  hemispherical  liners.  Extensive  use  has  been  made  of  other  liner 
designs  (discussed  later)  as  well  as  linear  and  circular  charges.  The  liner 
cross  sections  for  linear  and  circular  charges  are  wedges  and  semicircular 
configurations,  respectively.  The  jet  produced  by  a  linear  charge  is  in  the 
form  of  a  thin  ribbon  and  in  the  case  of  a  torus  (semi-circular)  charge  a 
tubular  spray  [28]. 

The  collapse  mechanism  of  various  shaped-charge  liners  ha3  been  verified 
via  flash  radiography.  These  novel  experimental  results  lead  to  the 
formulation  of  fundamental  Jet  collapse  theories.  These  theories  will  be 
discussed  later. 
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Figure  2.  Typical  Shaped-Charge  configuration. 


Figure  3.  Schematic  Collapse  of  a  Typical  Shaped-Charge 
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V.  APPLICATIONS 


Shaped-charges  are  extremely  useful  when  an  intense,  localized  force  is 
required  for  the  purpose  of  piercing  a  barrier.  The  main  application  is  in 
the  military  arena  including  torpedos,  missiles,  high  explosive  anti-tank 
(HEAT)  rounds  including  hand  held  (bazooka  type)  rounds,  gun  launched  (e.g., 
grenades)  cannon  launched  and  various  bombs.  The  targets  are  armors, 
bunkers,  concrete  or  geological  fortifications  and  vehicles.  Attacks  agains* 
aircraft  and  spacecraft  are  possible.  Underwater  applications  (torpedos)  are 
possible  with  the  design  such  that  water  does  not  enter  the  hollow  charge 
area.  In  fact,  most  warheads  of  the  type  described  above  contain  an  ogive  to 
cover  the  liner.  This  ogive  acts  as  an  aerodynamic  (or  hydrodynamic)  shield 
while  the  projectile  is  in  flight,  it  can  provide  a  housing  for  impact  fuzes 
or  guidance,  stability  and  control  electronics,  and  it  provides  a  built-in 
standoff  designed  to  improve  performance.  Figure  U  illustrates  the  concept  of 
standoff  distance  and  other  nomenclature  we  will  use.  The  standoff  is  the 
distance  from  the  base  of  the  charge  (or  liner)  to  the  target  in  question. 
Figure  5,  from  Mohaupt  [28]  shows  an  old  HEAT  artillery  projectile. 

The  largest  shaped— charge  application  was  known  as  the  Beethoven 
Apparatus  or  Histel  (mistletoe)  project  [2,533.  The  Mistel  technique  used  a 
fighter  aircraft  mounted  piggyback  on  the  top  of  a  large  bomber  aircraft.  The 
unmanned  bomber  carried  a  large  warhead  in  its  nose.  The  warhead  consisted  of 
a  2  m  (6  foot)  diameter,  wide  angle,  conical  shaped-charge.  The  liner 

probably  had  a  120°  apex  angle,  was  about  30  mm  thick,  and  made  of  either  mild 
steel  or  aluminum.  The  warhead  weighed  3500  kg  (7700  lbs)  and  the  explosive 
alone  weighed  1720  kg  (3800  lbs).  The  fighter  pilot  flew  the  combination  to 
the  target,  aimed  it,  released  it,  then  returned  to  his  base.  The  Germans 

developed  this  device  near  the  end  of  WW  II  and  most  of  them  were  captured 
intact. 

A  shaped-charge  with  an  aluminum  liner  was  used  in  an  attempt  to  place 
certain  identifiable  man-made  materials  into  orbit.  According  to  Kennedy  [2], 

a  35  included  angle  aluminum  shaped-charge  was  installed  on  a  multistage 
rocket  assembly.  The  craft  was  fired  into  near-space  and  the  warhead  was 
detonated  in  an  attempt  to  project  hypervelocity  fragments  into  earth  orbit. 
These  tests  were  conducted  at  Holloman  Air  Force  Base  in  1955-56.  It  was 
never  established  whether  or  not  any  aluminum  particles  were  detected. 

Along  the  same  line,  a  colleague  of  mine,  S.  K.  Golaski  of  the  Ballistic 
Research  Laboratory  and  the  old  Firestone  Tire  and  Rubber  Company  (now,  a  part 
of  Physics  International) ,  developed  shaped-charge  meteor  simulators  for  NASA 
15*4,55].  The  objectives  of  these  studies  were  to  obtain  the  luminous 
efficiency  of  a  meteor  like  body  of  known  mass,  composition  and  speed  during 
re-entry  into  the  earth's  atmosphere.  The  luminous  efficiency  is  the  percent 
of  kinetic  energy  of  the  body  which  is  converted  into  visible  light  as 
observed  by  photographic  study  of  the  visible  re-entry  tail.  The  study  was 
designed  to  obtain  the  required  mass,  composition  and  speed  from  pellets 
generated  by  a  specific  shaped-charge  design  flown  on  a  solid  state  propellant 
vehicle.  The  vehicle  was  used  to  achieve  the  necessary  velocities.  Nickel 
and  iron  were  the  shaped-charge  liner  materials.  The  rocket  engines  were 
designed  to  carry  the  shaped-charge  liner  above  the  atmosphere,  turn,  and 
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allow  the  detonated  shaped-charge  to  accelerate  at  hyper-velocity  through  the 
earth's  atmosphere  (re-entry).  The  intent  was  to  simulate  a  body  (meteor) 
re-entering  the  earths’  atmosphere.  A  specialized  shaped-charge  liner  design 
was  required  along  with  a  specialized  bi-explosive  waveshaping  device  [55]. 

Many  other  specialized  shaped-charge  applications  have  been  pursued  by 
the  Departments  of  Defense  of  several  nations.  These  specialized  designs 
included  confinement  or  tamping  of  the  explosive  fill,  varying  the  geometry 
and  material  used  for  the  confinement,  varying  the  geometry  or  type  of 
explosive  used,  altering  the  mode  of  initiation,  using  explosive  lenses  or 
more  than  one  type  of  explosive  or  an  explosive-non-explosive  barrier  or  gap, 
waveshaping  or  shaping  the  detonation  wave  (usually  done  to  insure  a  uniform 
wave  with  a  short  head  height,  see  Figure  H  for  definition  of  head  height),  or 
varying  the  standoff  distance.  Also,  significant  effects  can  be  achieved  by 
varying  the  liner  material  (Including  the  use  of  non-metals  such  as  glass), 
varying  the  liner  thickness,  varying  the  liner  diameter,  tapering  (or  causing 
a  gradual  wall  thickness  variation  either  continuously  or  discontinuously)  or 
varying  the  liner  geometry.  The  liner  geometry  variation  may  utilize  the  same 
basic  geometry,  e.g.,  varying  the  conical  apex  angle,  or  may  employ  a 
radically  different  liner  configuration.  Other  useful  liner  geometries  are 
hemispheres,  truncated  (from  the  equator)  hemispheres,  disc  or  dish  shaped 
(SFF  like)  devices,  tulips,  trumpets,  dual  angle  cones,  or  a  combination  of 
the  above  such  as  hemi-cones  or  tandem  devices.  In  fact,  any  arcuate  device 
may  be  used.  Also,  spin  compensated  liners  may  be  used,  especially  when 
associated  with  spinning  warhead  applications.  Spin  compensation  (i.e., 
causing  the  jet  to  spin  enough  to  compensate  for  the  spin  of  the  warhead)  may 
be  achieved  by  metallurgical  spin  compensation  or  by  the  use  of  fluted  liners 
(fluted  liners  contain  raised  ridges  either  on  the  outside  or  Inside  surface 
of  the  liner  to  allow  the  Jet  to  form  with  a  given  angular  momentum  to 
compensate  for  the  rotation  of  the  warhead  in  flight).  Without  spin 
compensation,  a  jet  formed  from  a  spinning,  detonated  warhead  will  exhibit 
radial  instabilities  (if  the  rotation  speed  is  high  enough)  and  dispense 
radially,  thus  reducing  the  effective  penetration  capability.  Spin  and  spin 
compensation  will  be  covered  in  more  detail  in  a  later  section. 

A  unique  shaped-charge  warhead  design  was  developed  by  Kennedy  and  others 
in  1967  to  1970  as  reported  by  Kennedy  [23.  This  particular  missile  design 
used  a  two-stage,  tandem  liner  designed  to  produce  a  precursor  or  pre-jet  to 
remove  the  guidance  and  control  package  located  in  the  ogive  of  the  missile. 
Otherwise,  the  main  Jet  would  have  to  penetrate  the  seeker  package  at  a  short 
(non-optimum)  standoff  distance.  Alignment  difficulties  in  the  tandem 
configuration  made  it  necessary  to  blend  the  precursor  liner,  its  standoff 
tube  and  the  main  conical  liner  into  a  single  piece.  The  result  was  the 
trumpet  design  configuration  which  is  used  today.  This  design  showed  improved 
penetration  with  a  relatively  small  standard  deviation  in  penetration. 

There  exists  numerous  other  tales  of  warhead  development,  the  point  being 
that  most  of  the  warhead  concepts  being  pursued  today  are  not  original 
concepts.  Other  warhead/liner  designs  will  be  discussed  later. 

Another  application  of  shaped-charges  is  in  demolition  work.  This  area 
has  both  military  and  industrial  application.  Buildings,  bridges  and  other 
structures  are  the  common  demolition  targets.  Usually  demolition  charges  are 
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hand  placed.  Demolition  charges  are  sometimes  constructed  and  placed  in  sites 
using  a  plastic  explosive  and  a  collection  of  liners.  This  technique  allows 
for  more  flexibility  and  adjustment  to  the  conditions  at  hand  than  a 
collection  of  fixed  charges.  Of  course,  experience  on  the  part  of  the  blaster 
is  required  and  this  technique  (the  adjustment  of  the  charge  to  the  task)  has 
delayed  bulk  charge  production  and  research  towards  this  particular 
application.  The  3haped-charge  principle  is  also  used  in  construction  work  to 
break,  crack  or  drill  holes  in  rock.  A  technique  known  as  raudcapping  is 
sometimes  used  to  break  rock  and  usually  utilizes  an  unlined  hollow  charge. 
Shaped— charges  are  also  used  in  construction  as  earth  movers,  in  tunneling,  or 
to  assist  in  well  drilling. 

Hollow  charges  are  also  used  as  a  source  of  earth  waves  for  geophysical 
prospecting  and  seismic  exploration.  Other  applications  occur  in  mining 
(surface  or  underground),  submarine  blasting,  timber  cutting,  breaking  log 
jam3,  breaking  ice  jams  and  steel  mill  furnace  tapping. 

A  Jet  tapper  is  a  shaped— charge  with  a  conical  metal  or  glass  liner.  The 
steel  mill  furnance  tapping  problem  is  depicted  in  Figure  6  [373.  The  tapping 
problem  requires  a  means  of  starting  the  flow  of  molten  steel  once  the  tap 
hole  has  been  plugged.  Steel  mill  jet  tapping  is  sometimes  called  salamander 
blasting  since  salamander  is  the  term  used  to  describe  the  large  mass  of  iron 
deposited  at  the  base  of  the  blast  furnance  after  it  has  been  in  operation  for 
some  time.  The  tap  hole  is  dug  out  as  far  as  possible  and  jet  tapper  is 
inserted  into  the  tap  hole  using  a  loading  pole.  The  jet  tapper  is  then 
detonated,  clearing  the  plug,  usually  with  a  high  rate  of  success  [373. 

A  commercial,  DuPont  type  jet  tapper  is  schematically  illustrated  in  Figure  7. 

The  shaped-charge  assembly  may  be  aligned  on  axis  or  offset  up  to  10  ,  as 
shown  in  Figure  7,  to  reduce  the  line-of-sight  thickness  of  the  plug  that  must 
be  perforated. 

Another  industrial  application  of  the  shaped-charge  is  in  the  internally 
coned  end  of  certain  detonators.  This  indented,  lined  cavity  acts  to 
concentrate  the  effect  along  the  axis,  recall  Bloera  [153. 

Cook  [563  used  the  original  Munroe  effect  to  engrave  or  stencil  letters 
onto  metal  plates. 

Other  examples  are  in  the  oil  well  industry  where  large  diameter,  but 
extremely  short,  lined  shaped— charges  are  used  to  penetrate  various  geological 
formations  to  increase  the  flow  of  oil.  Oil  well  perforation  problems  present 
extremely  difficult  design  problems  due  to  the  minimal  amount  of  allowable 
space  available  in  the  well,  the  short  standoff  distances  required  and  the 
hostile  environment  within  the  well  [573. 

Linear  (wedge  shaped,  V— shaped  or  W— shaped)  shaped-charges  are  also  used 
as  cutting  charges.  They  generate  a  ribbon  shaped  jet  used  to  cut  metals  and 
other  materials.  Commercial  cutting  charges  are  available  from  such  places  as 
Ensign-Bickford,  Explosive  Technology,  DuPont  and  others.  For  cutting  charge 
applications,  it  is  sometimes  advantageous  to  use  homemade  cutting  charges 
which  can  be  optimized  to  the  particular  problem  on  hand.  Cutting  charges  and 
hollow  charges  are  also  used  as  explosive  separation  devices,  as  bolt  cutters 
and  for  other  applications.  The  shaped-charge  effect  is  used  on  systems  for 
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separation,  deployment,  safety  destruct  and  similar  spacecraft  components. 
Garcia  [58],  Brown  [59]  and  Sewell  [60]  provide  additional  data  on  the  linear 
cutting  charge  concept.  Leidel  [61]  presents  a  discussion  of  the  design  of  a 
annular  cutting  charge. 

An  interesting  paper  by  Jones  [62],  discusses  the  perforation  of  artic 
sea-ice  by  large,  demolition  type,  shaped-charges.  The  following  table  was 
extracted  from  Reference  [62]  and  illustrates  the  charge  sizes  used  as  well  as 
the  penetration  depth  achieved  in  various  targets. 

General  references  on  application  are  found  in  Kennedy  [2],  Cook  [5],  The 
Blaster’s  Handbook  [37],  and  Bawn  and  Rotter  [631.  Certain  specific 
applications  of  interest  are  found  in  Cook  [56],  on  engraving;  Torrey  [64]  on 
a  review  of  shaped-charge  weapons  of  WW  II;  Huttl  [65]  on  underground 
blasting;  Clark  [66]  on  concrete  fragmentation;  Lawrence  [67],  and  Austin,  et 
al,  [68,69,70]  on  boulder  breaking,  drilling  and  seismic  exploration;  Draper, 
et  al,  [71]  on  drilling;  and  Davidson  and  Westwater  [72]  on  cratering  and 
boulder  breaking.  Excellent  overviews  are  available  in  Murphy  [733  and  Baum, 
et  al.  [74],  Extensive  studies  are  currently  underway  in  cratering  (earth 
removal),  boulder  breaking,  and  penetration  of  concrete  and  geological 
formations  at  several  Installations,  notably  the  Waterways  Experimental 
Station,  Vicksburg,  Mississippi.  The  demolition  and  cratering  work  performed 
at  the  Waterways  Experimental  Station  is  too  extensive  to  cover  here,  but 
their  studies  cover  many  of  the  applications  cited  earlier.  See  Joachim  [75] 
as  one  example. 

Additional  applications  relating  to  explosive  metal  interactions,  and 
which  require  an  understanding  of  the  Jet  formation  phenomena,  are  explosion 
welding,  explosion  calding,  or  explosion  forming  of  metal  parts.  The 
principles  of  explosive  welding  are  beyond  the  scope  of  this  text.  Basically, 
the  technique  involves  using  explosives  to  form  bimetallic  or  trimetallic 
parts.  A  tremendous  savings  in  material  usage  is  thus  possible  because  a  thin 
layer,  of  an  often  costly  material,  can  be  applied  to  a  cheaper,  load  bearing 
substrate.  The  thin  noble  and  expensive  material  is  usually  used  to  provide 
an  inert  or  corrosion  resistance  barrier  to  the  substrate  material. 
Applications  occur  in  chemical  vessels,  heat  exchanges,  desalination  plant 
piping  and  other  areas  where  caustic  environments  exist.  With  explosive 
welding  techniques  even  metallurgically  incompatible  metals  and  immiscible 
metal  combinations  can  be  bonded,  as  well  as  compatible  metal  systems. 
References  [76,  77  and  78]  provide  an  excellent  introduction  to  the  various 
aspects  of  explosion  welding.  Later  a  brief  discussion  of  this  important 
field  will  be  provided. 
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Figure  4.  The  Nomenclature. 
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Figure  5a.  High-Explosive  Antitank  Artillery  Projectile,  HEAT  {28]. 


V-7 


Figure  5b.  HFAT  (High  Explosive  Antitank)  Projectile. 


JET  TAPPER 


Figure  6.  Steel  Mill  Furnace  Tapping  [37]  . 
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Figure  7.  The  Jet  Tapper  Charge  (Offset)  [37]  . 


TABLE  I. 

FROM  J.  M  JONES  [62] 


DESIGN  AND  PERFORMANCE 

DETAILS  OF 

STANDARD  DEMOLITION 

SHAPED  CHARGES. 

No.  1,  6  in. 

No.  11,  30  lb. 

M2  A  3 

M3 

Total  weight  (lb.)  approx. 

10 

50 

15 

40 

Explosive  weight  (lb.) 

6.7 

27.5 

11.5 

30 

Charge  diameter  (in.) 

6.0 

12 

6.0 

9.5 

Charge  length  (in.) 

7.5 

11.75 

9.44 

15.5 

Design  Standoff  (in.) 

5.7 

19.0 

5.5 

15.0 

Casing  material 

steel 

steel 

fibre 

steel 

Cone  diameter  (in.) 

6.0 

12.0 

4.88 

9.0 

Cone  angie  (deg.) 

80 

60 

60 

60 

Cone  material 

steel 

steel 

glass 

steel 

Cone  thickness  (in.) 

.140 

.225 

.350 

.150 

Cone  weight  (lb.) 

— 

- 

1.7 

5.6 

Penetration 

Armour  depth  (in.) 

6 

- 

12 

20  + 

Average  diameter  (in.) 

— 

- 

1.5 

2.5 

Mild  steel  depth  (in.) 

9 

- 

- 

- 

Concrete  depth  (in.) 

30 

72 

30 

60 

Entrance  diameter  (in.) 

- 

6-8 

3.5 

5 

Minimum  diameter  (in.) 

2 

- 

2 

2 

Fresh  water  ice  depth  (in.) 

72 

144 

diameter  (in.) 

- 

- 

3.5 

6 

(standoff  42  in.) 

Permafrost  depth  (in.) 

72 

72 

(30  in. 

(50  in. 

standoff) 

standoff) 

diameter  (in.) 

- 

- 

6  -  1.5 

8-5 

Crater  diameter  at 

- 

- 

28 

28 

design  standoff  (in.) 
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VI.  The  Generalized  Equations 


The  general  conservation  equations  of  mass,  momentum  and  energy  may  be 
expressed  as 


ap 

at 


Du. 


"axT 


k: 

1  J 


> 


> 


and 


where 


3x. 

1 


(o. .u.) 
1 J  J 


ij 


6 . .  P,  or 
ij 


the  stress  tensor  equals  the  stress  deviator  tensor  plus  the  hydrostatic 
stress.  Also, 

Et  =  1/2  (uA  u.^)  +  E  or 

the  total  energy  equals  the  kinetic  energy  plus  the  internal  energy. 

The  substantial  derivative,  represents 


D_  _  _a_  a 

Dt  at  +  Ui  3x. 

i 

Thermal  energy  terms  do  not  appear  in  the  energy  equation  since  for  the 
micro-second  time  scale  involved,  there  is  no  time  to  radiate  or  conduct  heat. 

These  equations  constitute  the  basis  of  the  hydrocodes.  We  need  to 
relate  the  stress  tensor  to  hydrostatic  pressures  and/or  strains  and  strain 
rates.  Then  the  strains  and/or  strain  rates  must  be  related  to  the  velocities 
and/or  the  velocity  gradients.  The  equations  required  for  this  step  are  known 
as  the  constitutive  relationships  which  are  intended  to  convert  the  stress 
tensor  to  a  pressure  and  velocities  and  material  properties  (e.g.,  viscosity, 
shear  modulus,  etc.).  These  material  properties  are  hopefully  known,  or 
expressible  in  terms  of  known  quantities. 

Next,  it  is  necessary  to  specify  the  appropriate  initial  and  boundary 
conditions  for  the  partial  differential  equations.  It  is  also  necessary  to 
establish  explosive  burn  routines  and  provide  algorithms  to  couple  the  burn 
routine  to  the  liner,  for  example,  or  to  couple  two  complete  sets  of 
conservation  equations  such  as  for  a  penetrator  and  a  target.  Finally,  a 
numerical  solution  technique  and  a  stability  criterion  are  required.  This  is 
obviously  a  complex  procedure,  but  satisfactory  results  have  been  obtained 
from  several  available  hydrocodes.  Note  that  accurate  metal  and  explosive 
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equations  of  state  and  accurate  constitutive  equations  and  associated  material 
properties  are  prerequisite  to  a  good  solution.  Assuming  this  information  is 
available,  a  few  limiting  assumptions  will  be  made.  First,  the  following 
equations  are  available  for,  say,  the  behavior  of  a  metal  undergoing  explosive 
loading.  A  method  to  treat  the  explosive  is  assumed. 


Equation  Description 

Number  of  Equations 

Unknowns 

Continuity 

1 

4  -  P, 

u. 

1 

C\J 

II 

tI 

Momentum 

1  for  each  spatial  direction 

5  -  P, 

ui’  P 

i  =  1,  2, 

Energy 

1 

5  -  P, 

V  EI 

i  =  1,  2, 

Equation  of  State 

1 

3  -  P, 

P,  Ej 

Thus,  for  a  general,  three  dimensional  problem  there  are  six  equations: 

1  continuity,  3  momentum,  1  energy  and  1  equation  of  state.  There  are  also 
six  unknowns:  3  velocity  components,  the  density,  pressure  and  internal 


energy .  Various 

simplifying  assumptions  change 

the  situation  as 

shown 

below: 

Case 

Equations 

Unknowns 

3-D 

6 

6  ( p ,  P ,  Ej ,  u1 ) 

i  = 

2,  3 

2-D 

(axisymmetric ) 

5 

(one  less  momentum  equation) 

5  (p,  P,  Ej.,  ui) 

i  = 

1,  2 

Incompressible, 

2-D 

3 

(continuity  and  2  momentum 
equations;  p  =  constant) 

3  (P,  u±)  i  =  1, 

2 

Incompressible , 

3 

3  (P,  u±)  i  =  1, 

2 

2-D, 

steady-state 

(continuity  and  2  momentum 
equations;  p=  constant) 

No  initial  conditions 
required,  one  less 
independent  variable. 

The  axisymmetric  case  is  the  most  commonly 

used,  since  this 

is  usually  a 

valid  assumption  for  explosive-metal  interactions.  Note  that  a  tremendous 
simplification  is  possible  for  an  incompressible  flow  assumption.  However,  an 
incompressible  flow  implies  a  low  Mach  number  flow  and  the  complete  absence  of 
shock  physics.  This,  of  course,  is  not  realistic.  It  may,  though,  provide 
approximate  answers  as  will  be  seen. 
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Unknowns 


Equation  Description  Number  of  Equations 


Continuity 

1 

**  :  p  , 

u.  i  = 

l 

1,  2,  3 

Momentum 

1  for  each  spatial 
coordinate 

5 :  p , 

ui ,  P  i 

=  1,  2,  3  Add  P 

Energy 

1 

5  :  p, 

V  EI 

i  =  1,  2,  3  Add  Ej 

Equation  of  State 

1 

3  •  p » 

P,  Ej 

Close  System 

The  point  to  be  made  here  is  that  the  generalized  explosive-metal 
interaction  is  extremely  complex.  The  governing  equations  are  complicated; 
sufficient  boundary  and  initial  conditions  must  be  specified;  the  constitutive 
relationships  and  material  properties  under  the  intense  dynamic  conditions 
encountered  are  not  well  known;  equations  of  state  must  be  accurately 
specified;  and  numerical  algorithms  and  solution  techniques  must  be  accurately 
developed. 

Next,  further  simplications  will  be  presented  which  lead  to  analytical 
models  for  the  collapse  of  the  shaped-charge  liner.  The  next  section  starts 
with  the  one-dimensional  hydrodynamic  equation,  i.e.,  the  Bernoulli  equation. 


Case 

3-D 

2-D 

(axisymmetric) 


Equations 

6 

5 


Unknowns 

6  ( p  >  P »  Ej ,  u^ ) 

i  =  1,  2,  3 

5  (p,  P,  Ej,  u^)  Need  1  less  momentum 
i  =  1 ,  2  equation. 


Incompressible , 
2-D 


3  3  (P,  uA) 

Continuity  4  i  =  1 ,  2 

2  momentum. 


p  =  constant.  No  energy 
or  EOS  needed. 


Incompressible,  3  3 

2-D,  steady-state  As  Above  As  Above 


No  initial  conditions 
required,  one  less 
independent  variable. 
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DEFINITION  OF  TERMS 


-  spatial  coordinate 
t  -  time 

p  -  density 

-  velocity  components 
-  stress  tensor 

S,  .  -  stress  deviator 

1 J  ». 

ij  -  Kronecker  delta 

P  -  pressure 

Et  -  total  energy 

Ej  -  internal  energy 
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VII.  The  Bernoulli  Concept 


One  of  the  most  basic  equations  in  shaped-charge  jet  collapse  (and  in  jet 
penetration)  is  Bernoulli's  equation.  Due  to  the  repeated  application  of  this 
equation,  and  in  order  to  clearly  understand  the  assumptions  inherent  in 
Bernoulli’s  (and  modified  Bernoulli's)  equation,  it  will  be  derived  using 
several  methods.  The  intent  here  is  not  to  prove  "Theorems"  or  even  to 
introduce  rigor  in  our  methods,  but  to  establish  assumptions  inherent  to 
Bernoulli's  and  similar  equations. 

1 .  The  Energy  Equation 

Define  a  portion  of  fluid  between  sections  (1)  and  (2)  as  shown  below, 


u  being  the  velocity.  We  assume  one  dimensional  (1-D)  flow  and  follow 
standard  fluid  texts  such  as  Liepmann  and  Roshko  [ 1 ] . 

This  fluid  portion  can  be  considered  to  be  bounded  by  pistons  at 
(1)  and  (2)  instead  of  a  fluid.  These  pistons  are  equivalent  to  the  fluid 
they  replace,  and  help  to  clarify  the  work  done  on  the  system  which  consists 
of  a  definite  mass  of  fluid. 
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During  the  small  time  interval  in  which  the  fluid  is  displaced  to  the 
position  (1)'  -  (2)',  a  quantity  of  heat  q  may  be  added,  (as  external  heat). 

Tne  energy  law  yields, 

q  +  work  done  =  increase  in  energy, 

and  interpret  the  work  done  and  the  change  (increase)  in  energy  as; 

Work  done:  The  volume  displaced  at  (1)  is  the  specific  volume  (per  a  unit 
mass)  v^.  For  steady  flow  the  displacement  at  (2)  is  also  a  unit  mass,  or 

specific  volume,  v^.  The  work  done  by  the  pistons  during  this  displacement  is 

P  v  -  P  V  , 

1  1  2  2 

neglecting  any  external  work  between  (1)  and  (2). 

Change  in  energy:  The  flowing  fluid  has  internal  energy  e  and  kinetic 
2 

energy  1/2  u  ,  per  unit  mass.  After  the  displacement,  there  has  been  an 

2 

increase  of  energy,  e_  +  1/2  u„  corresponding  to  the  displacement  between  (2) 

d  2 

and  (2)'  and  a  decrease  of  energy,  e1  +  1/2  u^  from  (1)  to  (1)'. 

2  2 

Therefore,  increase  of  energy  =  (e2  +  1/2  )  -  (e^  +  1/2  u^  ) 

For  steady  flow,  the  energy  equation  becomes 

q  *  P1  V1  •  P2  V2  1  '2  '  81  *  ,/2  (U22  ‘  “l2) 
or  for  q  =  o,  the  adiabatic  energy  equation  yields 

P1  »,♦«,♦  1/2  u,2  •  P2  v2  ♦  -  1/2  u22  . 
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This  equation  relates  conditions  at  two  equilibrium  states  (1)  and  (2)  of  the 
flow.  They  are  valid  even  if  viscous  stresses,  heat  transfer,  or  any 
nonequilibrium  conditions  exist  between  (1)  and  (2).  However,  (1)  and  (2) 
themselves  must  be  equilibrium  states. 

If  equilibrium  exists  all  along  the  flow  (at  every  station),  the 
equilibrium  condition  Is  valid  continuously  and 

2 

P  v  +  e  +  1/2  u  =  Const. 

Further,  if  the  internal  energy  does  not  vary, 


Then 


2 

P  v  ♦  1/2  u  =  const . , 

and  if  the  density  is  constant  (v  =  — )  , 

2 

P  +  1/2  p  u  =  const . , 

which  is  Bernoulli's  equation  derived  from  thermodynamic  and  conservation  of 
energy  considerations. 

2.  Euler's  Equation 

Now  apply  Newton's  law  to  a  flowing  fluid,  where 
F  =  ma  . 

Now  view  the  flow  in  an  Eulerian  sense,  i.e.,  observe  the  flow  of  a  fluid 
particle  as  it  encounters  the  various  conditions  in  the  tube  through  which  it 
accelerates. 

The  acceleration  is  in  two  parts: 

Since  conditions  vary  along  the  tube,  there  is  a  velocity  gradient  in  the 

direction  of  the  flow,  The  rate  of  change  of  velocity  is  proportional  to 

this  gradient  and  the  speed  at  which  a  particle  moves  through  it.  Thus,  the 
acceleration  due  to  convection  through  the  velocity  gradient  is 

3u 

u  9x  ' 

Also,  conditions  at  a  given  station  may  be  changing  due  to  nonsteady 
(nonstationary)  effects  or 

au 

at  ’ 
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Then , 


9u  du 

a  =  —  +  u  — 
x  9t  ox 


Next,  it  is  necessary  to  determine  the  force  on  a  fluid  element  as  shown  below 


PA  +  aI(PA)  Ax 


A 


The  force  (Ax  A),  and  dividing  by  the  volume  of  the  element,  AAx,  and 

o  X 

dividing  through  by  the  density  p  gives  the  force  per  unit  mass  or 

=  _  1  3P 

p  3x 


This  result  is  valid  for  a  particle  (element)  or  any  shape, 
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converging,  diverging  or  whatever.  The  proof  is  standard  (e.g.,  [1])  using 
Gauss's  law.  Note  that  the  force  term  neglects  all  viscous  terms.  Then, 


3u  3u  1  dp 

at  u  3x  P  dx 

which  is  Euler's  equation. 

For  steady  flow,  u  ^  +  1  dp  =  0,  since  u  =  u  (x) ,  P  =  P  (x) ,  and  for 
1-D,  steady  flow.  dx  p 

For  incompressible  flow,  (p  =  p  )  and 
2 

P0  Y~  +  P  =  const. , 

which  is  the  classic  derivation  of  Bernoulli's  equation.  Now  consider  the 
momentum  equation  and  the  flow  through  a  control  surface,  where 


and 


u 


vpA)  +  u  —  (puA)  =  0.  (continuity  equation  x  u) 
Adding  these  equations  and  with  algebraic  manipulation,  we  have 

!_  (puA)  +  |_  (pu2A)  «  -  A  -  4t  (PA)  +  P  |A  , 


at  8x 

for  1-D  flow. 


3x 


3x 


Integrating  between  any  two  stations  (1)  and  (2),  yields 

2  2  2  2 
i  (PuA)  dx  =  p2u2  A2  -  p1u1  A2  =  PjAj  -  P2A2  +  P  dA  ,  and 


at 


2  pda  -  P^  (A2  -  A^),  defining  a  mean  pressure  P 


m. 
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The  integral  momentum  equation  is  very  general  in  that  it  is  valid  even  when 
there  are  frictional  forces  and  dissipation  regions  inside  the  control  volume, 
provided  they  are  absent  at  the  reference  stations  1  and  2.  The  integration 
of  the  differential  momentum  equation  corresponds  to  a  summation  of  the  forces 
on  adjacent  fluid  elements,  and  of  the  flow  into  them: 


The  forces  on  adjacent  internal  faces  are  equal  and  opposite,  and  thus  cancel. 
Similarly,  the  inflows  and  outflows  through  adjacent  faces  cancel,  leaving 
only  the  force  and  flux  at  the  control  space  boundaries.  If  there  is  a 
non-equilibrium  region  inside  the  control  space,  it  does  not  effect  the 
integral  result. 


For  steady  flow  in  a  constant  area  duct  (A^  =  A2), 


P2U2 


+  =  P. 


P1U1 


which,  sans  the  factor  of  1/2,  is  similar  to  Bernoulli's  equation. 

Harlow  and  Pracht  [2]  commented  on  the  momentum  conservation  and  the 
Bernoulli  principle  when  considering  jet  penetration.  They  observed  that  the 
conservation  of  momentum  principle,  assuming  the  jet  and  target  to  be  of  the 
same  width,  was  nearly  the  same  equation  as  derived  from  the  Bernoulli 
principle  where  the  target  is  effectively  infinite  in  width. 


3,  Euler's  Equation  Extended 


Euler's  equation  can  be  found  in  any  standard  text,  e.g.,  Shapiro  [3]. 
Euler's  Equations  of  motion  in  3-D  cartesian  coordinates  (as  derived  above  for 
a  1-D  case)  are  in  general: 


1  9P  9u  9u  9u  9u  Du 

p  9x  _  9t  +  U  9x  V  3y  +  W  9z  Dt 


O) 


I  !£. 

P  9y 


+ 


w 


9v 

9z 


D v 
Dt 


(2) 
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and 


I  il 

p  3z 


3w 

3t 


+  u 


3w 

9x  + 


Consider  2  cases: 


DW 

Dt 


(3) 


1.  Integration  along  a  streamline  {steady  flow)  [3].  A  streamline  is 
defined  by  the  condition  that  at  each  instant  the  velocity  vector  is  tangent 
to  the  streamline.  Thus, 


Now  multiply  equation  (1)  by  dx,  equation  (2)  by  dy  and  equation  (3)  by  dz, 

assume  steady  flow  (—  =  0)  and  for  a  particular  streamline  (using  equations 

(^),  3t 


3P 

3x 


dx  =  u 


3u 

3x 


dx  +  v 


3u  , 
t—  dx 

3y 


3u  , 

*  “  a?dx 


i 


from  equation  (1),  and  similiar  expressions  result  from  equations  (2)  and  (3). 
Then,  adding  the  resulting  equations  and  noting  that 


Jn  3P  ,  3P  ,  3P  , 
dP  =  87  dX  *  8?  dy  *  3?  d2> 
and 


,  3u  ,  3u  ,  3u  , 

du  *  87  dx  *  77  d>'  *  8?  dl-  etc 
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yields , 


-  dP 

P 


udu  +  vdv  +  wdw 


d 


2  2 
u  +  v 

2 


+ 


V  being  the  resultant  velocity. 
Thus , 


dP  =  -V  dV  and  for  incompressible  flow, 

V2 

p  +  P“  =  constant  along  a  streamline  =  Bernoulli  constant  of  the 
streamline. 

This  equation  is  only  valid  from  point  to  point  on  a  given  streamline.  Thus, 
the  Bernoulli  constant  may  vary  from  one  streamline  to  another.  Of  course,  if 
all  the  streamlines  originate  in  a  uniform  flow  region,  the  same  Bernoulli 
constant  would  prevail  over  the  entire  flow  field. 

2.  Integration  for  Steady,  Irrotational  flow  [3].  We  abandon  the 
restriction  that  only  points  along  the  same  streamline  are  considered  and 
assume  only  irrotational,  steady  flow.  This  means, 


5v 

3x 


3u 

3y 


m  .  sa  =  o 


3w 

3y 


dW  3v  ^  , 

^  °  ,  and 


3u 

3z 


3_w 

3x 


a u  ow  _  q 


Now,  multiple  equation  (1)  by  dx,  equation  (2)  by  dy,  and  equation  (3)  by  dz 
and  use  the  irrotationality  assumption,  which  yields 

1  3P  .  3u  ,  3u  ,  9u  .  ^  3v  ,  3w  d 

'  p  ^7  dx  =  u  dx  +  v  37  dx  +  w  3z  d  37  d  9x  d  d  ' 


etc.,  and  adding  the  equations, 


2  2  2  2 
—  f—  +  v  +  w  ,  _  3  X \ 

3x  2  3  "  3x  ' 


3 


etc.,  and  finally 


1  V 

-  i  dp  =  d  (y) 


which  is  identical  to  the  Bernoulli  equation  obtained  before.  Thus,  for 
steady,  continuous,  frictionless  and  irrotational  flow  the  Bernoulli  constant 
is  the  same  for  all  streamlines  and  relates  any  two  points  in  the  flow  field. 


The  constant  has  a  particular  value  for  each  streamline  except  for 
irrotational  flow  in  which  case  the  constant  is  the  same  for  all  streamlines. 
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4. 


Navier-Stokes  Equation 


We  begin  by  investigating  the  derivation  of  the  Navier-Stokes  equations 
following,  e.g.,  Schliohting  [4]. 


These  equations  represent  the  equations  of  motion  of  a  compressible, 
viscous,  Newtonian  fluid  in  three  dimensional  flow.  The  flow  field  will  be 

specified  by  the  velocity  vector  V  =  ui  +  vj  +  wk  and  the  equation  of  motion 
(Newton's  Second  law)  is 


=  F 
Dt 


+  P 


DV 

where  qj  is  the  substantial  or  Eulerian  derivative  of  the  velocity  vector  V. 

F  represents  the  body  forces  (forces  acting  through  the  mass  of  the  body) 
e.g.,  gravity,  and 

P  represents  the  surface  forces,  stresses,  shears,  friction,  and  pressure 
forces. 

We  ignore  F. 


The  surface  forces,  P,  depend  on  the  strain  rate  of  the  fluid.  The 
stress  is  related  to  the  rate  of  strain  via  an  isotropic,  Newtonian  fluid 
assumption.  This  means  that  the  relation  between  the  components  of  stress  and 
the  rate  of  strain  is  the  same  in  all  directions  (isotropic)  and  this 
relationship  is  linear  (Newtonian).  However,  shear  stresses  (viscous  effects) 
will  be  ignored. 


To  obtain  the  surface  forces  consider  the  fluid  volume  shown  below. 


Vol.  =  dx  dy  dz 
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3P 


The  x  (plane  to  x)  net  surface  force  component  is  dx  dy  dz  for  a 

one  dimensional  (1-D)  case.  The  resultant  surface  force  per  unit  volume  is 

3P  . 


3x 


where 


A 


Px  =  i  +  txy  j  +  tx2  k,  after  ijesolving  P  into  normal  stresses  (j-)  and 
shear  stresses  (t) . 


For  a  1-D  case  and  ignoring  the  shear, 

A 

P  =  i  — > 

«x 


and  neglecting  body  forces 

3o 


Du 

Dt 


3.x 


Hydrostatic  stress  systems  (fluid 
Stokes  hypothesis  will  not  be  invoked, 


pressure)  will  not  be  considered, 
but  deformations  and  rate  of 


deformation  will  be  neglected.  Thus,  the  normal  stress  will  not  be  related  to 
pressure,  velocities  or  velocity  gradients  as  done  in  the  classic 
Navier-Stokes  derivation. 


Instead  consider  a  uniform  stress 


field  acting  on  the  fluid  element  at 


rest  (assume  a  negative  o  )  somewhat  analogous  to  a  thermodynamic  pressure. 

Then,  expanding  the  substantial  derivative  and  for  incompressible  flow 

36 


,_3u  3u.  _  _ 

P  l3t  +  Sx'-  "  3x 


for  a  1-D  flow,  (recall  ^  =  ^+  u^j+v|£+wyi). 


Su1 


3u 


3u 


Dt 


at 


3xl 


Then,  for  a  1-D,  steady  state  system 
2 


d  ru  .  e  0 

p  a?  (rJ  *  dx  • 


or 


p  y-  +  a  =  const. , 


with  interpreted  as  a  uniform  normal 


3y 


3z 


stress.  This  result  is  analogous  to  a 
"Bernoulli  equation  with  strength"  which  with  the  basic  Bernoulli  equation  is 
frequently  cited  in  penetration  and  shaped-charge  jet  collapse  and  formation 
discussions.  Of  course,  the  full  Navier-Stokes  equations  relax  to  Euler's 
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equations  under  the  specified  assumptions,  which  in  turn  lead  to  the  classic 
Bernoulli  equation  as  shown  earlier. 

We  conclude  by  remarking  that  the  Bernoulli  equation  can  represent  a 
conservation  of  energy  or  a  conservation  of  momentum  depending  on  how  it  is 
derived.  We  also  note  that  the  "Bernoulli  with  strength"  type  equation,  or 
the  "modified  Bernoulli  equation"  has  some  theoretical  Justification. 
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VIII.  The  Gurney  Velocity  Approximation 

The  motion  of  a  metal  driven  by  an  explosive  was  studied  in  the  19^0's  to 
early  1950's  in  a  series  of  papers  by  Gurney,  Sterne  and  Thomas  [1,2, 3,  ^5]  as 
applied  to  the  motion  of  fragments  resulting  from  the  detonation  of  an 
adjacent  high  explosive  [6].  This  simple,  approximate  analysis  assumes  that 
the  potential  (or  chemical)  energy  of  the  explosive  charge  before  detonation 
is  converted  directly  to  the  kinetic  energy  of  the  metal  after  detonation  and 
to  the  expansion  of  the  explosion  products.  The  gas  detonation  products  are 
assumed  to  expand  uniformly  and  with  constant  density. 

The  so-called  Gurney  approximation  is  based  on  the  conservation  of 
momentum  and  energy.  The  results  represent  excellent  engineering 
approximations,  within  10},  of  the  experimental  results  or  detailed  numerical 
results.  The  Gurney  method  may  be  applied  to  any  one-dimensional 
explosive-metal  interaction  system.  The  amazing  accuracy  of  the  relatively 
simple  Gurney  formulae,  over  a  wide  range  of  metal  mass  (M)  to  explosive  mass 
(C)  ratios  (0.1  M/C  j<  10.0),  is  apparently  due  to  offsetting  errors.  These 
errors  are  caused  by  ignoring  rarefaction  waves  passing  through  the  detonation 
product  gases,  which  causes  the  calculated  velocity  to  be  too  high,  and 
assuming  an  initial  constant  density  distribution  of  the  detonation  product 
gases  rather  than  a  distribution  with  a  peak  at  the  surface  of  the  charge 
caused  by  the  detonation  wave,  which  causes  the  calculated  velocity  to  be  too 
low  [6].  Based  on  references  [6,7,8]  the  Gurney  formulae  will  be  derived. 

1 ♦  Open-Faced  Sandwich 

From  Kennedy  [8§],  a  specific  energy  (energy  per  unit  mass) 
characteristic  of  a  given  explosive  is  assumed  to  be  converted  from  chemical 
energy  in  the  initial  state  to  kinetic  energy  in  the  final  state.  The  final 
kinetic  energy  is  partitioned  between  the  driven  metal  and  the  detonation 
product  gases  by  an  assumed  linear  velocity  profile  in  the  gases. 

Figure  1  illustrates  the  open  faced  sandwich.  It  consists  of  a  slab  of 
explosive  confined  by  a  metal  on  one  side  only.  This  configuration  is  often 
used  in  experiments  to  obtain  constitutive  properties  of  materials  or  when  it 
is  necessary  to  impact  a  large  target  surface.  The  metal  plate  is  termed  the 
slapper  or  driver.  The  motion  of  the  metal  is  obtained  from  momentum  and 
energy  conservation.  For  the  open  faced  sandwich,  Figure  2  illustrates  the 
assumed  liner  velocity  distribution.  From  this  linear  distribution, 

VGAS  m  ■  <Vo  *  V>  r  -  V ■  (1) 

o 

The  y  coordinate  is  interpreted  to  be  Lagrangian  or  identified  with  a 
material  particle.  A  given  material  particle  is  assumed  to  move  at  constant 
speed  at  all  times  for  the  purpose  of  constructing  momentum  and  energy 
balances.  The  initial  thickness  and  density  of  the  explosive  are  y  and  p  , 
respectively.  The  energy  and  momentum  balances  follow:  °  e 

-vlV 


•Note  that  this  is  J.  E.  Kennedy,  not  the  Don  Kennedy  referenced  earlier. 
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Figure  1.  Open-Faced  Sandwich. 
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Figure  2.  Linear  Velocity  Distribution  for  the  Open-Faced 
Sandwich  Configuration  [8] . 
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where 


C 

E  =  Specific  explosive  energy 
M 
t 


Y  p 
o  e 


P  t 
m 


metal  thickness ,  and 


0  =  -  MV  +  p 


n 


(Vo  ♦  V) 

0 


dy  , 


and  finally 


V  =  i^2E 


M, 


t1*2?  M 

M„  +  C 


6(1  4  g) 


-1/2 


(3) 


(H> 


and 


Henry  [6]  provides  additional  mathematical  details.  Also,  Henry  presents  a 
complete  review  of  the  Gurney  formulae.  The  equations  for  some  common 
symmetric  and  asymmetric  configurations  are  presented  in  Figures  3,  ^  and  5. 

Figure  3  depicts  two  asymmetric  configurations,  namely  a  flat  sandwich  and 
a  cylinder.  M  denotes  the  metal  mass  and  C  denotes  the  explosive  mass.  These 
terms  are  obtained  as  the  product  of  the  density  and  thickness  of  the  material 
in  question.  Thus,  M  and  C  are  really  masses  per  unit  area.  The  term  V  is 
the  resultant  metal  velocity.  The  quantity  >/2E  occurs  in  all  Gurney 
equations.  It  has  the  units  of  velocity  and  is  termed  the  Gurney 
characteristic  velocity  for  a  given  explosive.  Gurney  characteristic 
velocities  have  been  tabulated  for  certain  explosives  by  Kennedy  [8],  Henry 
[6],  Jones  et  al  [7],  and  many  others,  notably  Dobratz  [9].  Several 
explosives  are  tabulated  from  Reference  [9]  with  their  density,  detonation 
velocity  and  Gurney  characteristic  velocity  or  Gurney  constant.  Note  that  two 
Gurney  constants  are  given,  one  for  warheads  in  which  confining  cases  rupture 
at  small  expansions  (prompt)  and  the  other  for  warheads  in  which  more  ductile 
case  materials  are  used  which  expand  further  before  rupturing  (terminal). 

The  tabulated  densities  are  the  actual  density  of  the  explosive  used  to 
determine  the  Gurney  energy.  The  detonation  velocities  sometimes  correspond 
to  a  different  explosive  density  than  that  given  in  the  table.  It  is 
recommended  that  the  explosive  of  interest  be  accurately  identified  and  the 
appropriate  parameters  be  taken  directly  from  Dobratz  and  the  references 
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Figure  3.  The  Flat  Sandwich  and  Cylindrical  Configurations. 


therein,  or  any  other  suitable  source  such  as  Meyer  [10]  or  Reference  [11]. 

La  Rocca  [12]  presents  a  simple  method  for  accurately  predicting  the  Gurney 
constant  and  the  relative  power  of  any  explosive. 

Figure  4  presents  the  spherically  symmetric  case  and  Figure  5  shows  two 
asymmetric  configurations,  an  open  faced  sandwich  and  an  asymmetric  sandwich 
which  has  two  metal  plates,  one  being  a  tamper.  The  derivation  of  these 
equations  is  relatively  straightforward  using  the  assumptions  stated  along 
with  elementary  calculus  and  algebra.  Derivations  are  given  in  Kennedy  [8], 
Jones,  et  al  [7]  and  Henry  [6].  Henry  [6]  considers  a  few  additional 
explosive/metal  geometries.  The  ratio  of  V//2E  or  the  ratio  of  the  metal 
velocity  to  the  Gurney  characteristic  velocity  is  an  explicit  function  of  M/C 
or  the  metal  to  charge  mass  ratio.  For  the  asymmetric  sandwich  where  the 
second  metal  plate  may  be  considered  to  be  a  tamper,  V/J2E  is  also  a  function 
of  N/C,  where  N/C  is  the  mass  to  charge  ratio  of  the  tamper  plate.  The 
velocity  expression  given  is  that  of  the  metal  plate,  M.  The  velocity  ratio 
is  plotted  as  a  function  of  M/C  in  Figure  6  for  all  the  non-tampered 
cases  presented.  Since  values  of  /2E  have  been  tabulated  for  certain 
explosives,  estimates  of  the  metal  velocity  can  be  obtained  for  a  given  M/C 
for  the  geometries  considered  here. 

Figure  7  plots  the  metal  velocity  as  a  function  of  M/C  for  an  open  faced 
sandwich  and  a  few  values  of  the  Gurney  constant,  /2E,  corresponding  to 
PBX-9404,  TNT  and  Comp  B.  This  plot  allows  easy  interpolation  or 
extrapolation  to  other  values  of  the  Gurney  constant. 

The  range  of  applicability  of  the  Gurney  formulae  is  restricted  due  to 
the  simplifying  assumptions  in  the  derivation.  These  restrictions,  according 
to  Kennedy  [8],  Jones,  et  al  [7]  and  Henry  [6]  are  listed  and  discussed  in  the 
following  tables.  The  Gurney  assumption  of  a  linear  velocity  profile  of 
constant  density  detonation  product  gases  greatly  deviates  from  conventional 
gas  dynamic  theory.  This  assumption  introduces  the  largest  error  in 
configurations  involving  a  free  explosive  surface  such  as  the  open  faced 
sandwich,  for  which  the  Gurney  method  may  overestimate  the  metal  velocity  and 
impulse.  Henry  [6]  used  a  parabolic  pressure  profile  in  the  detonation 
products  to  remove  the  assumption  that  the  gas  density  was  constant  at  any 
given  time.  A  simple  gas  equation  of  state  was  also  used.  Henry  concluded 
that  the  additional  mathematical  complexity,  inconsistent  with  the  basic 
simplicity  of  the  Gurney  method,  was  not  worthwhile  based  on  the  small 
improvement  obtained  in  accuracy.  Other  extensions  of  the  Gurney  method  are 
given  in  Reference  [13]. 
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Explosive 

Density 

D 

vfE",  km/s 

3 

gm/cm 

kra/s 

Prompt 

Terminal 

5  -  7  mm 

19  -  26  mm 

Comp  A- 3 

1.61 

8.47 

2.402 

1.59 

2.63 

Com  B 

1.71 

2.70 

1.717 

7.89 

2.35 

2.756  -  2.821 

1.717 

2.71 

Cast 

1.68 

2.402 

1.62 

2.32 

Pressed 

1.59 

2.335 

Comp  C-4 

1.52 

8.37 

2.176 

Comp  C-3 

1.60 

7.63 

2.68 

Cyclotol  75/25 

1.7514 

8.2  -  8.3 

Cast 

1.69 

2.286 

Pressed 

1.64 

2.362 

DATB 

1.68 

7.52 

1.975 

Explosive  D 

1.50 

6.85 

1.942 

HBX-1 

1.70 

7.31 

2.213 

* 
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Explosive 

Density 

,  .D 

v/2E,.  km/s 

gm/cm^ 

km/s 

Prompt 

Terminal 

5  -  7  mm 

19  -  26  mr 

HBX-3 

1 .81 

6.9  -  7.1 

1.984 

HMX 

1.89 

9.11 

2.97 

LX-14 

1.835 

8.83 

2.80 

NM 

1.14 

6.35 

2.41 

NQ 

1 .44 

7.65 

1.896 

Octol  75/25 

1 .81 

8.48 

2.58 

1 .821 

2.83 

DETASHEET  C 

(DuPont) 1 .48 

6.8 

2.1  -  2.3 

(BRL) 

7.0 

2.8 
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ExDlosive 

Density 

D 

»/5e,  km/s 

gm/cm^ 

kra/s 

Prompt 

Terminal 

5  -  7  mm 

19  -  26  m 

PBX-9011 

1.77 

8.5 

2.82 

PBX-9404 

1.84 

8.8 

2.90 

PBX-9502 

1 .885 

7.71 

2.377 

Pentolite  50/50 

Cast 

1.64 

7.52 

2.301 

Pressed 

1 .57 

2.317 

PETN 

1 .76 

8.26 

2.93 

RDX 

1.59 

8.25 

2.451 

1.77 

8.70 

2.93 

TACOT 

1 .61 

7.25 

2.12 

Tetryl 

1.63 

7.5 

2.274 

1 .62 

2.50 

TNT 

1.63 

2.039 

2.419  -  2.! 

1.63 

2.37 

Cast 

1.61 

6.73 

2.097 

Pressed 

1.54 

6.93 

2.103 
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Figure  4.  The  Spherical  Configuration. 


OPEN  FACED  SANDWICH 
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Figure  5.  Asymmetric  Configurations. 
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Figure  7.  V  Versus  M/C  for  an  Open-Faced  Sandwich,  Kennedy  [8]. 


RESTRICTIONS  ON  THE  GURNEY  MODEL 


Restriction 

Range  of  M/C  ratio 


Acceleration  Phase 


Direction  of  Detonation 
Propagation 


Gas  Velocity  Profile 
Assumption 


Remarks  and  Recommendations 


Henry  [6]  claims  accuracy  for  0.1  <  M/C  <  5  and 
Kennedy  [8]  states  a  range  of  0.2  <  M/C  <  10  for 
velocity  calculations.  Impulse  calculations  are 
usually  acceptable  for  M/C  >  0.2. 


The  Gurney  method  in  its  basic  form  is  not  capable 
of  analyzing  motion  during  acceleration.  However, 
Henry  [6]  and  Jones,  et  al  [7],  used  the 
assumption  of  a  liner  velocity  profile  in  the 
uniformly  dense  gases  in  conjunction  with  an 
equation  of  state  of  the  gases  to  yield 
acceleration  solutions  for  flying  metal  plates. 

The  detonation  product  gases  must  be  allowed  to 
expand  sufficiently  to  complete  the  acceleration. 
The  flying  metal  will  reach  its  calculated 
velocity  only  if  no  external  forces  (or 
interactions)  are  applied  during  the  initial 
acceleration  phase. 


The  detonation  of  the  HE  drives  the  metal  at  a 
given  velocity  (approximately)  for  a  given  M/C 
regardless  of  the  angle  between  the  detonation 
front  and  the  metal.  The  direction  in  which  the 
metal  is  driven  will  vary  slightly  with  the  angle 
(see  Taylor  angle  discussion). 


The  assumed  linear  velocity  profile  and  constant 
density  gas  expansion  are  major  assumptions. 
Ignoring  effects  of  both  rarefaction  waves  and  of 
pressure  peaks  near  the  metal  surface  would  seem 
to  be  cancelling  errors.  However,  these 
assumptions  are  inherent  to  the  simplicity  of  the 
Gurney  method. 
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RESTRICTIONS  ON  THE  GURNEY  MODEL 


Restriction 

One-Dimensional  Motion 

Metal  Strength  Effects 

Metal  Spallations 

Early  Case  Fracture 


Remarks  and  Recommendations 


The  Gurney  approach  is  not  valid  for  the  estimate 
of  variations  in  local  velocity  of  a  plate  driven 
by  a  charge  with  a  tapered  thickness  (or  a  tapered 
metal  thickness).  In  this  case,  an  average  value 
of  M/C  could  be  used  to  estimate  the  final 
velocity  of  the  entire  plate. 


Forces  exerted  by  the  metal  to  oppose  deformation 
are  not  considered,  other  than  inertia.  Hoop 
stresses  in  cylinders  and  spheres  can  reduce  the 
metal  velocity  in  explosions  and  the  strength 
effect  is  greater  in  implosions. 


Metal  spallation  may  occur  when  M/C  <  2  for  high 
density  explosives  and  metals.  Spallation  can,  in 
come  cases,  be  avoided  by  introducing  an  air  gap 
of  a  few  millimeters  between  the  explosive  and  the 
metal.  This  results  in  a  small  decrease  in  metal 
velocity. 


The  leakage  of  the  detonation  product  gases 
through  fractures  in  the  metal  case  can  decrease 
the  final  metal  velocity  perhaps  significantly 
depending  on  the  nature  and  extent  of  the  leakage. 
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The  specific  impulse  delivered  by  the  explosive  can  be  calculated  as  the 
total  momentum  imparted  to  the  metal  body  divided  by  the  total  explosive 
charge  mass  or 


I 

sp 


Mv 

C 


(6) 


The  specific  impulse,  I  ,  will  now  be  derived  for  an  unconfined  surface 

charge.  Explosive  may  be  detonated  directly  on  the  surface  of, a  massive  metal 
body  in  order  to  deliver  a  desired  impulse  for  testing  pruposes.  If  the 
loaded  body  is  considered  to  be  rigid,  all  detonation  product  gases  will  flow 
away  from  the  surface  and  a  maximum  specific  impulse  will  be  delivered  to  the 
body.  For  a  very  large  metal  body  mass  to  explosive  mass  ratio,  M/C  >>  1,  the 
open  faced  sandwich  velocity,  equation  (4),  reduces  to 


V 


then  I  % 
sp 


3E 

2 


Also  since 


C 

M 


,  from  equation  (6). 


V  |^2  ^  +  lj  from  equation  (5)  and  for  M/C  »  1, 

2V  M/C  and  with  V  given  by  equation  (7), 

2  /2E  or 

/6E 


(7) 


(8) 


for  the  maximum  free  gas  velocity. 

The  impulse  delivered  to  a  body  by  the  explosive  can  be  increased  by 
tamping  or  confining  the  explosive  with  an  outer  layer  of  metal.  Tamping  acts 
to  hold  back  the  detonation  gases,  keep  the  pressures  high  and  force  the 
detonation  gases  to  do  useful  work  against  the  body.  A  tamping  and  a 
confinement  are  analogous.  The  Gurney  method  can  be  used  to  estimate  the 
change  in  impulse  as  a  function  of  the  tamper  areal  density.  If  the  heavy 
metal  body  sandwich  is  considered  to  be  rigid,  it  resembles  the  plane  of 
symmetry  in  the  symmetrical  flat  sandwich  configuration  shown  in  Figure  3. 

The  velocity  of  the  tamper  plate  is  then  given  by  the  flat  sandwich  of  Figure 
3  with  the  tamper  to  explosive  mass  ratio  N/C,  used  instead  of  M/C  or 


V  +  /2E 


(9) 
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The  effective  specific  impulse  of  the  explosive  for  a  tamping  ratio  of  N/C 
follows  from  the  equation  for  the  asymmetric  sandwich  of  Figure  5,  namely, 


V 

m 


/2E 


1  +  A" 


3(1  +  A) 


.2  M 
A 


-1/2 


(10) 


Assuming  M/C  >>  1  and  with  A  = 


1  +  2  M/C 
1  +  2  N/C 


,  the  specific  impulse  becomes 


1/2 

I  ,  tamped  =  /2E  (N/C  +  1/2)  /  (N/C  +  1/3) 
sp 


(11) 


Note  that  when  N/C  =  0,  i.e.,  no  tamper,  equation  (11)  reduces  to 


I  as  derived  earlier, 

sp  12 

The  acceleration  of  explosively  driven  metal  plates  can  also  be 
determined  by  using  the  Gurney  assumption  of  a  liner  velocity  profile  in 
uniformly  dense  gases  in  conjunction  with  an  equation  of  state  for  the  gases. 
Acceleration  calculations  showing  the  displacement  as  a  function  of  time  or 
speed  are  important  when  estimating  the  effects  of  venting  in  a  perforated 
casing.  Henry  [6]  and  Jones,  et  al  [7],  calculate  the  acceleration  of  the 
casing  for  various  sandwich  configurations.  These  same  authors  also  extended 
the  Gurney  analysis  by  considering  the  pressure  gradients  within  the  gas  in 
order  to  remove  the  assumption  that  the  gas  density  is  constant  at  any  given 
time.  These  results  however,  complicate  the  simple  Gurney  model. 

The  Gurney  model  yields  explicit  algebraic  relationships  for  estimating 
the  velocity  imparted  to  a  metal  in  contact  with  a  detonating  explosive.  The 
model  can  be  used  for  several  simple  geometric  metal/explosive  configurations. 
The  specific  impulse  of  an  explosive  can  also  be  obtained  for  these  geometries 
and  is  directly  related  to  the  Gurney  energy  of  the  explosive.  Thus,  the 
Gurney  method  can  be  directly  and  simply  applied  to  a  multitude  of 
explosive/metal  interaction  problems.  Design  and  parametric  studies  can  be 
readily  performed.  Applications  include  warhead  and  fragmentation  design,  the 
study  of  the  efficiency  of  conversion  of  the  high  explosive  chemical  energy  to 
the  kinetic  energy  of  the  plate,  explosive  initiation  by  the  impact  of  an 
explosively  driven  plate,  the  launching  of  a  dielectric  plate  by  electrical 
detonation  of  a  metal  foil,  and  the  calculation  of  layers  of  metal  fragments 
in  conjunction  with  shock  wave  physics. 

Note  that  the  Gurney  method  does  not  consider  the  propagation  of  shock 
waves  through  metals,  but  only  the  terminal  effect  of  the  shock  propagation. 
Duvall  [14]  provides  an  excellent  description  of  shock  wave  propagation, 
transmission  and  reflection.  He  describes  a  plane  wave  lens,  a  plane  wave 
generator  and  methods  of  forming  oblique  shock  waves.  He  also  states  the 
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pressure  levels  attained  by  the  various  methods  of  shock  wave  generation. 
Duvall1 s  article  [14]  explains  shock  wave  phenomena  without  recourse  to 
advanced  mathematical  concepts.  Thus,  Reference  [14]  is  an  excellent 
introductory  article  for  the  layman.  Courant  and  Friedrichs  [15]  provide  an 
excellent  mathematical  treatment  of  shock  wave  physics.  Shock  wave 
propagation  through  solid  media  is  an  important  aspect  of  explosive/metal 
interaction,  but  is  beyond  the  scope  of  this  introductory  report. 

An  extension  of  the  Gurney  method  was  given  recently  by  Chanteret  [ 1 6 ] • 
Chanteret  developed  an  analytical  model  for  symmetrical  geometries  with  energy 
partitioning  during  expansion.  This  model  allows  a  calculation  of  the  Gurney 

energy  as  a  function  of  the  expansion  ratio.  The  expansion  ratio  is  (R/R  )n 

o 

where  R  is  the  current  value  of  the  radius  of  the  explosive-metal  interface 
Rq  is  the  undetonated  radius  of  the  explosive-metal  interface.  The 

constant  n  depends  on  the  geometric  configuration  (n  =  1  for  symmetrical  plane 
sandwiches  and  n  =  2  for  cylinders  with  a  constant  wall  thickness). 

Chanteret' s  results  can  be  used  in  conjunction  with  the  classical  Gurney 
formulae  to  obtain  the  entire  velocity-time  curve  for  all  symmetrical 
steady-state  configurations.  Also  Chanteret  extended  the  Gurney  model  to 
include  imploding  cylinders  by  introducing  a  fictitous  rigid  boundary  in  the 
explosive  cylinder.  The  results  are  claimed  to  agree  well  with  available 
experimental  data  and  two-dimensional  hydrocode  results. 

The  Gurney  method  and  extensions  of  this  method  are  useful  analytical 
tools.  The  flyer  plate  velocity  depends  only  on  the  M/C  ratio  and  the  Gurney 
constant  or  Gurney  characteristic  velocity,  /2IT.  This  constant  is  available 
for  several  explosives  (as  we  have  seen),  but  when  the  Gurney  constant  is  not 
available,  Kennedy  [8]  recommends  that  E  “  0.7  be  used  when  is  the  heat 

of  detonation.  This  recommendation  follows  from  the  fact  that  0.61  <  EHD  < 

0.76  for  most  explosives.  Another  option  is  to  use  La  Rocca's  method  [12]  to 
calculate  the  Gurney  constant. 

The  Gurney  method  ignores  rarefaction  effects  and  the  shock 
characteristics  of  metals.  Thus,  we  can  expect  errors  to  occur  when  these 
effects  are  important.  An  example  would  be  a  short  open  ended  cylinder  which 
will  deviate  from  the  Gurney  prediction  near  its  ends.  Also,  the  difference 
in  shock  impedance  between  different  metals  will  alter  the  energy  coupling  and 
as  a  result,  the  same  M/C  ratios  with  different  metals  but  with  the  same 
explosive  will  yield  different  velocities.  The  applicable  range  of  M/C  values 
is  0.2  <  M/C  <  10  as  discussed  earlier.  Note  that  as  M/C  decreases,  e.g., 
when  adding  HE  to  a  given  metal  slab,  the  velocity  cannot  increase 
indefinitely  as  predicted  by  the  Gurney  method. 

Finally,  the  Gurney  equations  assume  that  the  metal  moves  in  a  direction 
normal  to  its  surface.  This  is  true  when  the  detonation  wave  encounters  the 
metal  at  normal  incidence.  It  is  not  true  when  the  detonation  wave  encounters 
the  metal  plate  at  a  grazing  incidence.  The  grazing  incidence  is  treated  by 
another  model  known  as  the  Taylor  angle  approximation. 

Before  investigating  the  Taylor  angle  approximation,  a  few  Gurney  type 
formulae  relating  to  flat  plates  and  imploding  devices  will  be  discussed. 
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For  shaped  charges,  a  popular  formula  for  a  single  flat  plate  backed  by 
slab  of  explosive  is  given  by  [17],  where  P  &  M/C,  and  V  is  the  liner 
collapse  velocity,  and 


V' 


V2 


3  u 


1/2 


(12) 


For  this  same  configuration,  with  the  detonation  wave  (of  velocity  D) 
propagating  in  a  direction  tangent  to  the  liner,  Duvall  and  Erkman  [18,  19] 
used  hydrodynamic  theory  to  obtain  the  collapse  velocity  formula 


Kleinhanss  [20]  also  presented  equation  (13)  and  credited  Trinks  with  another 
formula  based  on  semi-empirical  data  related  to  the  explosive  launching  of 
plates , 


2  , 

V  =  0.36D  arctan  C“=T7>  *  ( 1  ^ ) 

o 

These  formulae  were  developed  for  flat  plates,  and  although  they  are 
quite  accurate,  they  do  not  account  for  the  effect  of  curvature  such  as  in 
implosive  geometries,  e.g.,  a  conical  shaped  charge.  Kleinhanss  [20],  in 
experiments  with  imploding  cylindrical  charges,  showed  a  strong  dependence  of 
the  liner  collapse  velocity  (VQ)  on  its  radius.  He  obtained  the 
semi-empirical  formula 


_ 1 _ 

[CQ  +  e  f (b) ) 


(15) 


where  e  is  the  metal  thickness,  b  *  rQ  -  r^,  rQ  and  r^  being  the  outer  and 

inner  explosive  radii,  respectively,  and  thus  b  is  the  explosive  thickness. 

Co  and  f(b)  are  empirical  constants  that  depend  on  the  type  of  explosive  and 

metal  being  used.  Other  studies  on  implosive  geometries  were  given  by 
Chanteret  [16]  as  discussed  earlier,  Chou  et  al,  [21]  and  Hennequin  [22].  All 
of  these  studies  provide  improved  results  over  the  basic  Gurney  model  for 
implosive  geometries. 

Other  formula  are  available  for  plate  acceleration  and  are  often  used  in 
shaped  charge  liner  collapse  and  explosive  welding  models.  For  example,  Singh 
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[23]  considers  the  velocity  of  the  casing  of  a  thin  walled  cylinder  filled 
with  explosive  and  detonated  at  one  end.  The  velocity  of  the  casing  at  the 
time  of  fracture  is  given  by 


=  I  (1  ♦  sinO 
D  ■+ 

where  each  element  of  the  casing  subtends  an  angle  a  with  the  axis  of 
propagation  of  detonation. 

Deribas  [24]  gives 


(16) 


w  ,  l 


jin 


x  i  +  44-  +  l 

^  t 

where  n  is  the  charge  to  mass  ratio. 
Shushko,  et  al  [25],  gives 


(17) 


TT=  °-61 


t\~  +  I2n  *  IS 


+  ltn  *  9 


n(n  +  6) 


(18) 


Mikhailov  and  Dremin  [26]  provide  several  frequently  used  expressions  for 
the  plate  speed.  They  are: 


V 

o 


0 . 6  02- 
2  +  n 


V 

o 

D 


S  -  1  _  Hff 
c  "7  Dt 


(19) 


(20) 


(21) 


where  y  = 


16  1 
27  n  ’ 


0  = 


t  is  the  time  taken  to 


travel  the  flight  distance  or  the  length  of  the  driver  plate,  and  H  is  the 
height  of  the  explosive  charge. 
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Smith,  et  al  [27]  used  formulae  analogous  to  equations  (12)  and  (17)  to 
determine  the  terminal  velocity  of  flyer  plates  used  in  plate  acceleration 
experiments.  Smith,  et  al  [27]  also  calculated  the  flyer  plate  terminal  bend 
angle  which  will  be  discussed  shortly.  Both  the  bend  angles  and  flyer  plate 
terminal  velocities  are  in  agreement  with  the  experimental  values  obtained  for 
three  different  explosives. 

Explosive-metal  interaction  models  are  also  useful  in  fragmentation 
calculations.  Fragmentation,  of  course,  does  not  always  result  from  an 
explosive-metal  interaction.  For  example,  fragments  are  generated  when  gas 
pressurized  pipes  or  vessels  fail.  The  velocity  of  the  fragments  following  a 
gas  pressurized  system  failure  is  necessary  in  order  to  determine  the  damage 
to  neighboring  equipment  on  structures.  Fragment  velocities  resulting  from  a 
pressurized  pipe  or  vessel  failure  are  usually  predicted  by  the  Moore  equation 
[28].  The  Moore  equation  is  semi-empirical  based  on  experimental  data  from 
the  velocity  of  explosively  generated  fragments  (i.  e.,  Gurney  data). 

Another  method  for  estimating  the  fragmentation  velocity  from  ruptured 
pipes  or  vessels  follows  from  the  theoretical  model  of  Baum  [29], 

Of  course,  other  Gurney  like  models  are  also  available.  The  section  on 
explosive  welding  and  the  sections  on  jet  formation  provide  references  to 
other  metal  acceleration  models.  Walters  [30]  and  Walters  and  Harrison  [31] 
provide  several  references  (mostly  Soviet),  related  to  the  acceleration  of 
metals  by  high  explosives.  Dehn  [32],  recently  published  a  report  containing 
an  overview  of  Gurney  type  velocity  approximation  models  as  well  as  the  Taylor 
angle  (grazing  incidence)  models.  Dehn  provides  additional  historical 
information  and  correctly  analyzes  geometric  configurations  in  this  report 
such  as  the  Jelly  Roll,  the  Dagwood  and  similar  explosive-metal  multilayer 
arrangements. 

In  conclusion,  care  must  be  taken  in  comparing  the  various 
explosive-metal  interaction  formulae.  Many  authors  use  M/C  (ratio  of  the  mass 
of  metal  to  the  explosive  mass)  as  a  key  parameter.  Other  authors  prefer  C/M 
as  the  parameter. 

Next,  grazing  detonation  waves  will  be  investigated. 


VIII-21 


THE  TAYLOR  ANGLE  APPROXIMATION 


The  Taylor  approximation  [33]  is  shown  in  Figure  8,  which  depicts  a 
grazing  (parallel)  incidence  of  the  detonation  wave  to  the  metal  surface.  The 
plate  is  deflected  at  an  angle  0  from  its  initial  position.  Acceleration  of 
the  metal  to  its  final  velocity  is  assumed  to  be  instantaneous.  For  this 
steady-state  condition,  the  metal  plate  is  assumed  to  undergo  pure  rotation, 
i.  e.,  no  net  shear  flow,  or  no  change  in  length  or  thickness.  Thus,  the 
plate  element  that  was  at  P  initially  will  be  a  P'  after  launch  and  the 
lengths 


OP  =  OP' 


Construct  a  line  from  0  perpendicular  to  PP'.  This  line  bisects  the  angle 
since  OPP*  is  an  isosceles  triangle.  If  the  time,  t,  is  measured  from  the 
time  when  the  detonation  wave  passes  point  P,  then: 

OP  =  Dt , 


PP'  =  Vt 


and  .  0  pp'/2 

sm  -r  =  — - 

OP 


Vt 

2Dt 


3L 

2D 


(1) 


An  estimate  of  the  direction  in  which  the  metal  plate  is  projected  is  the 
angle  0/2  which  can  be  determined  from  equation  (1)  with  V  known  from  the 
Gurney  method  and  D  determined  from  the  specified  explosive.  As  the  plate 
moves,  it  will  be  tilted  an  an  angle  0  from  its  Initial  position,  and  by 
assumption,  will  not  be  rotating. 


When  using  a  streak  camera  located  perpendicular  to  the  initial  charge 
axis,  a  component  of  the  velocity,  namely  Va,  the  apparent  velocity,  was 
measured  by  Kury,  et  al  [31*].  This  velocity  can  be  related  to  V  by  the 
geometry  or 


Va  =  D  tan  6  , 


then 


V  _  2D  sin  6/2 
Va  D  tan  0 


2  sin  6/2  cos  9  (2  sin  6/2) 


sin  8 
cos  6 


2  sin  2  cos 


e 

2 


(2) 


(3) 


V  _  cos  8 
or  Va  ~  cos  0/2 


Hoskin,  et  al  [13]  expressed  their  results  in  terms  of  the  velocity 

component  Vvl,  the  velocity  component  normal  to  the  flight  plane  of  the  plate, 
N 
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Figure  8.  Direction  of  Metal  Projection  By  a  Grazing  Detonation  Wave. 


since  this  velocity  is  of  interest  when  the  flying  plate  is  used  in  impact 
studies.  Since  the  angle  between  VA  and  is  0,  we  can  write 


=  D  sin  6,  and 


(*) 


V_  =  2D  sin  6/2 

VM  D  sin  6 
N 


cos  0/2 


sec  0/2 


(5) 


The  difference  between  V,  and  is  usually  only  a  few  percent.  Also,  the 

term  V/2D  is  approximately  the  same  for  many  explosives.  Thus,  0  is 
approximately  constant.  If  we  fix  the  mass  of  the  plate  (its  thickness)  being 
driven,  we  observe  that  increasing  the  detonation  velocity  of  the  explosive 
implies  an  increase  in  the  plate  velocity  such  that  V/2D  is  approximately 
constant.  These  concepts  will  be  used  to  study  jet  formation  theory. 
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IX.  JET  FORMATION 


Birkhoff,  Mac  Dougall,  Pugh  and  Taylor  [1]  formulated  the  first  theory  of 
conical  shaped  charge  jet  formation  assuming  that  the  detonation  wave  produces 
such  large  pressures  during  the  collapse  process  that  the  material  strength  of 
the  liner  may  be  neglected.  In  fact,  the  liner  is  treated  as  an  inviscid, 
incompressible  fluid.  The  conical  liner  is  modeled  as  a  wedge,  and  a  steady 
state  collapse  model  is  assumed.  Thus,  the  liner  elements  are  instantaneously 
accelerated  to  their  final  collapse  velocity.  The  steady  state  model 
predicted  a  Jet  with  a  constant  length  equal  to  the  slant  height  of  the  cone. 

It  has  been  observed  however,  that  shaped  charge  jets  possess  a  velocity 
gradient  with  the  tip  traveling  much  faster  than  the  tail,  causing  the  jet  to 
stretch  and  eventually  break  up.  This  steady-state  theory  was  later  modified 
by  Pugh,  Eichelberger  and  Rostoker  [2]  to  include  the  jet  velocity  gradient. 

The  modified  "non-steady"  theory  is  based  on  the  same  principles  as  the 
original  theory  except  that  the  velocities  at  which  the  various  liner  elements 
collapse  is  not  constant,  but  depends  on  the  original  position  of  the  element 
in  the  liner.  We  begin  with  an  examination  of  the  Birkhoff,  et  al  [1]  steady 
state  model. 

As  the  detonation  wave  impacts  the  liner  the  pressure  on  all  sides  of 
the  liner  is  assumed  to  be  equal  and  the  liner  walls  are  assumed  to  collapse 
inward  at  a  constant  velocity,  say  VQ.  The  angle  26  between  the  moving  walls 

is  greater  than  the  original  apex  angle  2a,  because  of  the  finite  time 
required  for  the  detonation  wave  to  sweep  the  liner  surface  from  apex  to  base. 
Figure  1  describes  the  geometry  of  the  collapse  process.  The  conical  liner  is 
symmetric  so  a  represents  one  half  of  the  liner  apex  angle  and  6  represents 
the  collapse  angle.  The  liner  assumes  a  velocity,  VQ,  which  bisects  the  angle 

APP'  of  Figure  1.  To  show  this  consider  a  coordinate  system  having  a  constant 
velocity  such  that  the  origin  moves  from  point  P  to  P'  in  a  unit  time.  In 
this  coordinate  system  steady  state  conditions  exist  at  the  origin  and  the 
liner  moves  inward  along  P’P  and  flows  outward  along  the  path  PA.  The 
velocity  of  the  liner  passing  through  this  region  changes  its  direction  but 
not  its  magnitude,  because  the  pressure  forces  are  everywhere  perpendicular  to 
the  motion.  P'P  and  P'B  represent  the  entering  and  emerging  velocities, 
respectively,  of  the  liner  in  the  moving  coordinate  system.  P'B  is 
constructed  parallel  to  PA  and  P'P  and  P'B  are  equal  in  magnitude.  Since  the 
velocity  in  the  moving  coordinate  system  is  PP'  the  velocity  of  the  collapsing 
liner  in  the  stationary  system  is  the  vector  sum 

PP'  +  P’B  =  PB  =  V  . 

0 

Since  P'P  =  P'B,  the  triangle  BPP'  is  isosceles  and  since  P'B  is  parallel 
to  PA,  angle  BPP'  =  angle  PBP'  =  angle  BPA .  Call  these  angles  6.  Thus,  VQ 
bisects  angle  APP'. 

The  walls  of  the  collapsing  liner  are  two  planes  moving  inward.  The 
junction  of  these  planes  moves  from  A  to  B  with  a  velocity  V^, 

=  Vq  cos  (— 2“! - )  /sin  6  . 
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From  Figure  1,  the  angle  OPA  =  <P  follows  from  the  fact  that  a  +  $  +  180  -  g  = 

1 80  or  4>  =  8  -  i.e.f  the  sum  of  the  angles  of  a  triangle  equals  180°. 

Next,  consider  the  triangle  PXB  where  180  =  90  +  6  +  ($XBP  =  y)  or  y  =  90  -  6. 
In  triangle  XAB,  1 80  =  90  +  8  +  }ABX  or  $ABX  «  90  -  B.  Then  along  the  line 


OPP',  26  +  $  =  180  or  e  =  90  -  <}>/2  =  90 


Y  +  6  =  90,  y  +  90 


90  or  y  -  o  |  Angle  PP'B  = 

triangle  PP*B  we  have  26  +  iPP’B  =  180  or  1 80  -  <J>  +  }PP'B  =  1'80  or  $PP’B 

This  establishes  the  geometry  of  Figure  1.  Now  for  some  trigonometry. 


■m 


,  from  above.  Then  since 
d>  since  in 


In  triangle  APB  by  the  law  of  sines,  V./sin  6  =  V  /sin  8  or  V.  =  V  sin  6/ 

\  I  o  10 

sin  B.  The  sin  0  =  sin  (90  ^  ~~  j  =  003  ^6  -  a  ^  and 

V-,  =  V  cos  ( —  rr  -)  /sin  B  =  V  cos  (4>/2)/sin  B  .  ^ 

^  0  \  L  '  O 


A  moving  observer  positioned  at  A  would  observe  any  point  P  in  the  upper 
plane  approaching  him  with  a  velocity  cos  3  +  Vq  cos  9.  This  velocity,  V2, 

becomes  V,  =  V.  cos  B  +  V  cos  (90  -  =——))=  V.  cos  g  +  V  sin  (-  ~  — ) ,  and 

from  equation  (2),  \  '  2  /  "  ; 


V  =  V 
1  0 


cos 

. 


/ tan  B  +  sin 


(3) 


Also  Up  =  U  cos  a  where  U  is  along  PP'  and  using  the  law  of  sines  for  triangle 


PBP  • , 


V 


u 


sin  (6  -  a)  ”  sin  9 
V 


or 


U  = 


sin  (8  -  a) 


sin  90 


■M- 


V  cos 
o 


(S-2 


sin  (B  -  a) 


then 


Un  V  cos 
D  0 


M 


cos  a  sin(B  -  a) 


(4) 
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Returning  to  V2  from  equation  (3)  the  observer  sees  a  "jet"  moving  to  his 

right  and  a  "slug"  moving  to  his  left,  as  in  Figure  2.  Also,  as  viewed  by 
this  observer,  the  whole  process  appears  to  be  unchanged  by  the  lapse  of  time, 
or  steady  state  motion  exists.  This,  with  the  earlier  assumptions  of 
inviscid,  incompressible  and  one-dimensional  flow  allows  the  use  of  the 
Bernoulli  equation  or 

2 

p  =  Jj  1)  =  constant,  (5) 

which  relates  the  pressure  and  the  corresponding  velocity  U.  The  pressure  at 
any  point  in  the  liner  determines  the  velocity  at  that  point.  Assume  the 
liner  moves  away  from  the  detonation  so  fast  that  the  pressure  of  the  surface 
decreases  rapidly  and  the  pressure  on  all  surfaces  of  the  collapsing  liner  is 
constant.  This  allows  free  streamlines,  i.e.,  the  boundary  streamlines  are  at 
constant  pressure  and  density  (and  hence  velocity).  As  viewed  by  the 
observer,  the  jet  and  slug  will  appear  to  recede  with  exactly  the  same  speed, 
V2,  as  the  walls  approach  as  in  Figure  1. 

Returning  to  the  stationary  system  of  coordinates,  it  is  seen  that  the 
jet,  traveling  to  the  right  in  Figure  2,  has  a  velocity 

V  =  V  +  V  , 

while  the  slug,  traveling  to  the  left  in  the  moving  system  of  Figure. 2, 
has  a  velocity  to  the  right  given  by 


To  further  visualize  this  process  consider  that  the  point  P  (fixed  in  the 
upper  plane)  travels  to  point  B  (fixed  in  space)  in  unit  time,  the  material 
from  the  inner  surface  of  the  upper  plane  included  between  PA  and  AB  moves 
into  the  jet,  and  the  front  of  the  Jet  moves  to  the  right  a  distance  PA  +  AB 
in  the  same  time,  forming  a  high  velocity  Jet.  This  jet  velocity  is 
V  =  V1  +  V2  and  from  equations  (2)  and  (3): 


Also,  the  outer  surface  of  each  plane  forms  a  slug  moving  with  the  lower 
velocity,  Vfl  =  V1  -  V2  or 
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Observer  Stationed  at  the  Moving  Juncti 


Conservation  of  momentum  dictates  the  division  of  material  between  the  jet  and 
the  slug.  Let  m  be  the  total  liner  mass  per  unit  length  of  the  two  collapsing 
planes  (sides  of  the  wedge)  approaching  the  junction.  Let  m^  represent  that 

part  of  the  liner  mass  entering  the  jet  and  ms  represent  that  part  entering 
the  slug,  so 


m  =  m , 


m 


(8) 


Equating  the  horizontal  momentum  components  entering  and  leaving  the 
junction  A  of  Figure  2  in  the  moving  coordinate  system  yields 


"  V2  cos  8  ’  "s  v2  -  V2  .  (9) 

Solving  equations  (8)  and  (9)  simultaneously  gives 

mj  =  (m/2)  U  -  cos  6)  ,  ..  (10) 

and  * 

ms  =  (m/2)  (1  +  cos  6)  (11) 

According  to  this  model,  the  velocities  of  the  Jet  and  slug  and  their 
cross-sectional  areas  are  constant. 

The  analysis  assumed  a  wedge  configuration.  A  conical  liner  may  be 
treated  in  the  same  way.  In  the  conical  case,  the  walls  converge  on  the  axis 
from  all  sides.  The  moving  observer  must  travel  at  the  same  rate  as  in  the 
case  of  the  wedge.  In  order  for  the  process  to  appear  stationary  to  him,  the 
total  mass  per  unit  distance  along  the  axis  must  be  constant.  This  is 
approximately  true  for  a  constant  wall  thickness  liner.  In  reality,  the  liner 
wall  thickness  would  have  to  be  tapered  to  be  inversely  proportional  to  the 
distance  from  the  apex. 


For  the  detonation  wave  traveling  parallel  to  the  axis  with  speed  U^,  as 
in  Figure  1,  equations  (6)  and  (7)  become,  from  equation  (4), 


V 


cos  a 


sin  (3  -  a) 


esc  6  +  ctn  S  +  tan 


3  -  a 


(12) 


and 


V, 


cos  a 


sin  (3  -  a) 


esc  3  +  ctn  3  +  tan 


(13) 


The  Jet  velocity,  V,  increases  as  a  decreases,  since  3  also  decreases.  The 
jet  velocity  approaches  a  maximum  as  a  -*■  0,  or  V  =  Up  {1  +  cos  +  tan  8/2}  or 


B  +  0  as  a  +  0  so  V  =  2  U^,  which  means  that  the  jet  velocity  cannot  exceed 
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twice  the  detonation  velocity.  Also,  as  a  8  ^O,  Vg^O.  Note  that  as 

a  0  the  conical  liner  approaches  a  cylinder.  Cylindrical  lines  are  capable 
of  generating  high  velocity  (and  low  mass)  jets. 

For  the  hypothetical  case  of  a  conical  wave  front  moving  perpendicular  to 
the  surface  of  a  conical  liner  such  that  the  wave  strikes  all  surfaces  at  the 
same  time,  6  =  a  and  the  velocities  of  the  jet  and  slug  from  equations  (6)  and 
(7)  this  time,  become: 


V  = 


sin  a 


(1  +  cos  a),  and 


V, 


V 

o 

sin  a 


(1  -  cos  a)  . 


With  this  type  of  wave  front  the  jet  velocity  could  be  increased  indefinitely 
by  decreasing  a.  However,  as  “  ->  o,  VQ  ■+  0,  m^  “*•  0  and  the  jet  momentum 
mV 

nij  V  =  — ^ —  sin  a  -*•  0. 

Finally,  the  theory  of  Birkhoff,  et  al  predicts  steady  state  jet  and  slug 
velocities  as  in  equations  (6)  and  (7)  or  equations  (12)  and  (13)  and  masses 
as  in  equations  (10)  and  (11)  for  conical  or  wedge  shaped  liners. 

Comments  on  the  Birkhoff,  et  al  Solution 

The  steady  state  model  provides  quantitative  agreement  with  flash 
radiograph  experiments.  The  model  tends  to  overpredict  the  jet  velocities. 
Also,  the  jet  possesses  a  velocity  gradient  from  tip  to  tail  and  continues  to 
stretch  after  the  walls  have  collapsed  giving  a  jet  length  greater  than  the 
slant  height  of  the  cone.  The  ductility  of  the  jet  under  the  intense 
pressures  and  dynamic  collapse  conditions  is  responsible,  at  least  in  part, 
for  the  elongation  of  the  jet. 

Important  modifications  have  been  made  to  the  steady  state  theory  of 
Birkhoff,  et  al  notably  by  Pugh,  Eichelberger  and  Rostoker  [2]  and  Godunov, 
Deribas  and  Mali  [3].  Pugh,  et  al,  hereafter  referred  to  as  the  PER  theory, 
developed  a  non-steady  theory  based  on  the  same  concepts  as  Birkhoff,  et  al 
[1]  except  that  the  collapse  velocities  of  the  various  liner  elements  are  not 
the  same  for  all  elements  but  vary  depending  on  the  original  position  of  the 
liner  element.  Thus,  the  collapse  velocity  decreases  continuously  from  the 
cone  apex  to  its  base,  producing  significant  jet  elongation. 

In  the  Soviet  Union  Godunov,  et  al  [3]  hereafter  referred  to  as  the 
visco-plastic  theory,  modified  the  steady  state  theory  to  include 
visco-plastic  behavior,  i.e.,  strain  rate  dependence. 

Both  of  these  models  will  be  discussed  in  turn. 
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The  PER  Theory 


The  PER  theory  assumes  a  variable  instead  of  a  constant  collapse  velocity 
for  the  walls  of  the  conical  (or  wedge)  liner.  This  condition  greatly 
improves  the  steady  state  theory  results.  The  collapse  velocity  decreases 
from  the  apex  to  the  base  of  the  cone  and  Figure  3  illustrates  the  effect  of 
these  velocity  variations.  Note  that  as  the  collapse  angle  3  increases,  the 
jet  velocity  decreases,  but  the  portion  of  the  liner  entering  the  jet 
increases.  Figure  3  also  shows  how  the  decreasing  collapse  velocity  increases 
6.  As  the  detonation  wave  travels  from  P  to  Q  along  the  conical  surface  APQ, 
the  liner  element  originally  at  P  collapses  to  J.  The  liner  element 
originally  at  P',  starts  later  and  collapses  slower  than  P,  and  arrives  at  M 
at  the  same  time  P  arrives  at  J.  The  element  P’  would  have  reached  N  when  P 
reached  J  if  their  collapse  velocities  were  identical.  Thus,  with  a  constant 
collapse  velocity,  the  surfaces  remains  conical,  or  QNJ  is  a  straight  line. 
However,  since  P'  has  a  slower  collapse  velocity  than  P,  the  collapsing  liner 
has  the  non-conical  contour  QMJ  of  Figure  3*  The  angle  6  is  greater  than  B+ , 
the  "steady-state"  3.  This  assumes  that  each  liner  element  is  thin  and  not 
affected  by  its  neighbors,  consistent  with  the  hydrodynamic  assumption. 

Next,  Figure  4,  which  is  similar  to  the  steady  state  model  of  Birkhoff, 
et  al,  illustrates  the  flow  situation.  QJ  is  parallel  to  PA  and  equal  in 
length  to  PQ.  If  the  magnitudes  of  QP  and  QJ  are  equal  to  U,  they  represent 
the  velocities  in  the  moving  coordinate  system  of  the  liner  elements  entering 

and  leaving  the  region  P.  The  vector  PJ  =  VQ  is  the  velocity  of  the 

collapsing  liner  element  in  a  stationary  system  of  coordinates.  Thus,  the 
arguments  of  Birkhoff,  et  al  [1]  are  valid.  The  liner  element  does  not  move 
perpendicular  to  its  original  surface,  but  along  a  line  that  makes  a  small 
angle  6  with  the  normal.  This  angle,  from  Figure  4,  is 

■  x  v  n\\  (14) 

sin  6  =  Vo/ZU 

If  V  is  constant,  6  =  B+  and  =  (ft  — )  and  all  the  PER  results  relax  to 
the  steady  state  model. 

By  the  proper  choice  of  coordinate  systems,  the  geometric  relationships 
at  the  moving  junction  (J)  are  shown  in  Figure  5.  The  axis  of  the  cone  is 
along  JR  and  OJ  is  the  element  of  the  liner  moving  toward  the  axis.  This 

element  has  a  velocity  OJ  =  V  in  moving  coordinates.  The  velocity  of  the 

moving  coordinate  is  JR  s  Using  the  law  of  sines  from  Figure  5  we  have 

V  =  V  cos  (a  +  6) /sin  3,  and  (15) 

o 

V,  =  V  cos  (3  -  a  -  6)/sin  £. 

1  0  (16) 
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Figure  3.  The  Collapse  Process  for  a  Variable  Collapse  Velocity  Liner, 
PER  [2]. 
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^STAGNATION  POINT' 


Figure  5.  Relationship  Between  VQ,  The  Liner  Collapse  Velocity,  vjhe 
Collapse  Velocity  Relative  To  The  Collision  Point,  and  Vj  The 
Collision  Point  Velocity. 
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In  fixed  coordinates,  the  velocities  of  the  jet  and  slug  are  V  =  +  V  and 

V  =  V,  -  V,  respectively.  Using  equations  (15)  and  (16),  along  with  some 
si 

trigonometric  manipulations,  the  jet  and  slug  velocities  become: 


V.  =  V  (esc  6/2)  cos  (a  +  6  -  6/2)  and 
3  ° 

=  Vq  (sec  6/2)  sin  (a  +  6  -  6/2)  ,  ^g^ 

respectively.  Note  that  for  6  =  8+  =  a  +  26  we  get  the  Birkhoff ,  et  al  Jet 
and  slug  velocities,  after  some  trigonometric  manipulation. 

We  can  eliminate  from  equations  (17)  and  (18)  using  equation  (14)  to 

yield 

(19) 

Vj  =  Vo  ^CSC  6/2^  cos  -  8/ 2  +  sin-1  Vq/2u)  , 

and  Vg  =  V  (sec  6/2)  sin  (a  -  8/2  +  sin'1  V  /2u)  . 

o 

(20) 


These  equations  are  valid  in  either  the  steady-state  case  when  VQ  is  constant 

or  in  the  nonsteady  case  where  VQ  varies.  In  the  steady  state  case,  however, 

8  could  be  expressed  in  terms  of  a »  U  and  VQ  and  need  not  appear  in  the 
equation. 

Of  the  four  unknowns  V.,  V  ,  m.  and  m  we  still  need  expressions  for  m 

J  s  j  s  j 

and  m  .  These  equations  follow  from  the  conservation  of  mass  and  momentum, 

or,  where  m  is  the  liner  mass  and  m,  and  mg  are  the  jet  and  slug  masses 
respectively,  J 


dm  =  dm 


j 


+  dm  , 
s 


.  2 
sin 


(6/2),  and 


(21) 


dms 

dm 


cos2  (6/2)  . 


(22) 


Equations  (21)  and  (22)  are  identical  to  the  steady  state  theory,  (equations 
(10)  and  (11)). 

Equations  (19),  (20),  (21)  and  (22)  describe  the  mass  separation  and  the 
velocities  of  each  element  of  a  conical  liner.  They  depend  upon  the  cone 
angle  2a  and  the  detonation  velocity  Up  =  U  cos  a,  and  on  the  collapse  angle  8 

and  velocity  V0>  both  of  which  are  variables. 
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From  Figure  3,  observe  that  the  calculation  of  the  angle  g  +  would  be  more 
straight  forward  than  the  calculation  of  g ,  since  VQ  is  different  for  each 

element  of  the  liner  and  thus  the  contour  of  the  collapsing  liner  is  not 
straight  but  curved  as  line  QMJ  of  Figure  3.  Toward  this  end,  let  the 
cylindrical  coordinates  of  M  be  (r,  Z)  and  the  coordinates  of  P1  be  (X  tan  a, 
X)  in  Figure  3.  The  coordinate  directions  are  indicated  in  Figure  6.  Then  we 
have 


X  ♦  V  (t  -  T)  sin  A,  and 

(23) 

X  tan  a  -  VQ  (t  -  T)  cos  A  , 

(24) 

where  t  is  the  elapsed  time  since  the  detonation  wave  passed  the  apex  of  the 

cone,  T  =  =  X/(u  cos  a)  and  the  angle  A  =  a  +  6.  The  slope  of  the  contour 

UD 

of  the  collapsing  liner  at  any  time  t  can  be  obtained  from  the  3r/SZ  using 
equations  (23)  and  (24)  and  some  calculus.  The  time  when  a  given  element 
reaches  the  axis  is  found  from  equation  (24)  with  r  =  0  or 

„  T  _  X  tan  a  , 

V  cos  A  (25) 

o 

and  9  r/  3Z  evaluated  at  r  =  0  is  the  tan  g.  Using  the  Taylor  relationship , 
equation  (14),  to  simplify  we  have  finally 

sin  a  +  2  sin  6  cos  A  -  X  sin  a  (1  -  tan  A  tan  61  v'/V 
tan  g  =  - — _ _ _ _ _ o  o  .  (26) 

cos  a  -  2  sin  6  sin  A  +  X  sin  a  (tan  A  +  tan  6)  v'/V 

0  o 

Since  26  =  g+  -  a  and  2A  =  g+  +  a>  equation  (26)  may  be  further 
simplified  to  yield, 


sin  g+  -  X  sin  a  (1  -  tan  A  tan  6)  v'/V 
tan  g  =  — - _  o  o  .  (27) 

+  t 

cos  g  +  X  sin  a  (tan  A  +  tan  6)  V  /V 

o  o 

The  term  V  represents  the  partial  derivative  of  V  with  respect  to  X. 

0  o 

+  t 

The  angle  g  >  g  since  V  <  0,  i.e.,  the  collapse  velocity  decreases 

from  apex  to  base,  and  for  cone  angles  ( 2a )  which  are  not  extremely  large. 
Additional  mathematical  details  and  a  more  complete  discussion  can  be  found  in 
Reference  [2]  and  Eichelberger 's  Doctoral  dissertation  [4]. 

The  reader  should  pay  careful  attention  to  the  quadrant  in  which  the 
angles  calculated  in  equations  (26)  and  (27)  may  fall.  For  angles  in  the 


IX-13 


second  quadrant,  minor  modifications  of  the  equation  may  be  necessary  and 
numerical  algorithms  that  calculate  trigonometric  functions  should  be  checked 
for  the  proper  quadrant. 

The  PER  theory  provides  expressions  for  m^,  ms,  Vj,  Vfi,  6  and  E,  although 

the  number  of  unknowns  exceed  the  number  of  independent  equations. 

Experimental  evidence  that  verifies  the  model  and  the  assumptions  of  the  PER 
theory  was  given  by  Eichelberger  and  Pugh  [5].  Also,  Eichelberger  [6]  noted 
a  systematic  discrepency  between  the  theoretical  and  experimental  results.  It 
was  noted  that  the  liner  acceleration  to  the  axis  was  assumed  instantaneous  in 
the  PER  model  but  the  effect  of  a  finite  acceleration  time  would  influence  the 
value  of  6.  In  1962,  Allison  and  Vitali  17],  in  another  classic  paper,  used  a 
radioactive  tracer  technique  and  obtained  excellent  agreement  between  the  PER 
theory  and  experiment. 

Although  Allison  and  Vitali  obtained  excellent  agreement  between  the 

theory  and  experiment  for  a  42°  unconfined  copper  cone,  they  were  unable  to 

get  the  PER  theory  to  agree  with  the  data  from  a  42°  confined  copper  cone. 
Merendino  and  Jonas  of  the  Ballistic  Research  Laboratory  obtained  similar  data 

for  a  60°  unconfined  copper  cone  which  also  failed  to  match  the  results  of  the 
PER  model.  G.  H.  Jonas  recalls  that  Allison  concluded  that  the  PER  model 
provides  valid  predictions  for  shallow  angle,  unconfined  conical  charges,  but 
needs  to  be  modified  for  wide  angle  or  confined  concial  charges. 

The  PER  theory  constitutes  the  basis  of  all  one-dimensional  shaped  charge 
jet  formation  models.  Many  extensions  have  been  made  to  the  PER  model  and 
auxiliary  equations  have  been  added  to  facilitate  and  extend  the  solution. 

The  plate  bending  angle  or  the  Taylor  angle  ( 6)  can  be  determined  by  the 
formula  proposed  by  Richter  [8]  and  Defournaux  [9]: 

—  -  —  +  P  £  K 

26  e  *  (28) 

where  p  and  e  are  the  density  and  thickness  of  the  liner  wall,  respectively, 
and  e  is  the  explosive  thickness  driving  the  liner,  K  and  $  are  constants 

which  are  determined  from  the  type  of  explosive  used  and  the  angle  of 
incidence,  i,  which  the  detonation  wave  makes  with  the  liner.  Equation  (28) 
is  coupled  with  the  PER  theory  in  the  one-dimensional  Jet  formation  codes  BASC 
[10,  11)  and  DESC  [12]  and  probably  others.  Defournaux  [9]  also  presents  a 
good  discussion  of  the  liner  collapse  and  formation  which  is  an  excellent 
supplement  to  the  PER  study.  Kerdraon  [13,  14]  modified  equation  (28)  by 
including  terms  to  account  for  the  effect  of  liner  radius  in  imploding 
cylindrical  liner  charges. 

The  PER  theory  assumes  the  liner  element  is  instantaneously  accelerated 
to  the  axis  as  depicted  in  Figure  7(a).  The  first  order  correction  is  the 
assumption  that  the  liner  acceleration  is  constant  over  a  finite  period  of 
time,  as  first  proposed  by  Eichelberger  [6]  and  used  in  Reference  [12].  The 
velocity  then  increases  linearly  over  a  short  period,  as  shown  in  Figure  7(b), 
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until  it  reaches  its  final  velocity  or  collapses  on  the  axis.  In  Reference 
[12],  the  acceleration  is  given  by 


where  PcJ  is  the  Chapman-Jouguet  pressure  of  the  explosive,  p  and  c  are  the 

liner  density  and  thickness,  respectively,  and  c  is  an  empirical  constant. 
Figure  7  was  taken  from  Chou  and  Flis  [15]. 

Possibly,  a  more  realistic  form  of  the  liner  velocity  history  is  given  a3 
an  exponential  form  by  Randers-Pehrson  [16], 


as  shown  in  Figure  7(c).  This  formula  requires  knowledge  of  the  time 
constant,  t,  for  which  Chou,  et  al  [17]  recommend 


where  M  is  the  initial  mass  per  unit  area  of  the  liner  and  c.  and  c?  are 
empirical  constants.  1 


Considering  the  effects  mentioned  above,  Chou,  et  al  [17,  18]  provide  a 
more  exact  expression  for  the  projection  (Taylor)angle, 


1  •  1  i 

7TV0  Vc 


where  the  prime  denotes  differentiation  along  the  liner.  Randers-Pehrson  [ 1 6 ] 
proposed  a  similiar  formula.  Both  References  [ 1 6 ,  17,  18]  show  better 
agreement  with  the  hydrocode  results  than  the  Taylor  angle  (steady  state) 
formula  for  6. 


Models  of  this  type,  for  acceleration  to  the  liner  axis  during  the 
collapse  process,  have  been  incorporated  into  the  PER  theory.  These  one 
dimensional  computer  programs  use  the  analyses  presented  in  this  section  to 

calculate  the  Taylor  angle  6,  the  flow  velocity  V,  the  collapse  velocity  V  , 

»  ° 

the  stagnation  point  velocity  ,  the  collapse  velocity  gradients  V  ,  the  Jet 

o 

velocity  the  slug  velocity  V  ,  the  collapse  angle  8,  the  Jet  mass  m.,  and 
J  “  J 

the  slug  mass  m  .  The  input  values  consist  of  the  complete  geometry  and 
s 

material  density  of  the  conical  liner  and  the  explosive  charge.  The  mode  of 
initiation  of  the  charge  and  the  type  of  HE  (high  explosive)  must  be  known  to 
determine  the  HE  density  and  detonation  velocity.  The  mode  of  initiation 
determines  the  shape  of  the  detonation  wave  (i.e.,  plane  or  spherical)  and  the 
angle  of  incidence  of  the  detonation  wave  to  the  liner  wall.  Also,  the 
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relative  distribution  of  the  kinetic  energy  entering  the  jet  and  slug  are 


cos2  (a  +  5  -  6/2)  (29) 


sin2  (5  +  6  -  /2)  (30) 


respectively. 

The  stretching  jet  radius  is  also  calculated  assuming  the  Jet  to  be 
circular  in  cross-section.  In  addition,  the  material  near  the  liner  apex  does 
not  have  sufficient  time  to  reach  its  theoretical  collapse  velocity- and  hence 
the  earliest  formed  portions  of  the  Jet  have  lower  velocities  than  the  jet 
material  formed  behind  them.  This  "piling  up"  of  jet  particles  causes  the 
formation  of  the  inverse  velocity  gradient  and  this  coalescence  of  particles 
is  treated  separately  to  calculate  the  Jet  tip  velocity.  All  these  effects 
are  included  in  most  jet  formation  codes  such  as  BASC  [10,  11]  along  with 
certain  other  parameters  such  as  momentums,  jet  strains  and  strain  rates. 

To  calculate  the  strain,  jet  length  and  jet  radius,  it  is  necessary  to 
relate  the  position  of  a  Jet  element  back  to  its  original  position  in  the 
liner.  This  was  done  by  Carleone  and  Chou  by  using  a  simple  Lagrangian 
coordinate  definition  [19].  The  coordinate  x  is  used  for  the  axial  liner 
position  and  the  coordinate  £  for  the  position  of  the  jet  as  shown  in  Figure 
8.  The  Jet  position  is  given  as 


given  as 

dE. 

dE 

and 


dE 

dE~ 


C  (x,  t)  =  Z  (x)  +  (t  -  tQ)  V.  (x),  {31) 

where  t  is  the  time  when  the  liner  element  first  arrives  at  the  axis,  Z  (x) 
o 

is  the  location  of  the  formation,  and  Vj  (x)  is  the  Jet  velocity.  Note  that 

Vj  is  a  function  of  a  liner  position  but  £  is  a  function  of  both  liner 

position  and  time.  Outside  the  inverse  velocity  gradient  region,  a  liner 
element  with  a  smaller  x  always  proceeds  (is  ahead  of,  or  has  a  larger  £  than) 
an  element  with  a  larger  x.  The  strain  rate,  v,  is  given  as 


n 


Lim 


AC  -  CQ  (H/3x)t 

AC  =  (3C/9X). 

to 


1 


(32) 


Equations  (31)  and  (32),  with  the  statements  made  earlier,  allow  the  jet 
length  and  radius  to  be  calculated. 

Special  techniques  are  necessary  to  calculate  the  jet  tip  velocity  as 
discussed  earlier.  The  PER  theory  indicates  that  the  jet  velocity  decreases 
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Figure  8.  Relations  Between  the  Liner  Coordinate  x,  and  the  Jet  Coordinate  f 
Chou  and  FI  is  [15]. 


monotonically  from  tip  to  tail.  Thus,  the  gradient  of  the  Jet  velocity  is 
negative  with  respect  to  the  original  liner  position  or  positive  with  respect 
to  the  jet  position.  As  mentioned  earlier,  liner  elements  near  the  apex 
collide  on  the  axis  before  reaching  their  final  collapse  velocity,  VQ.  This 

reduced  collapse  velocity  results  in  a  reduced  jet  speed,  so  that  the  jet 
velocity  gradient  with  respect  to  the  liner  position  is  positive.  This  is 
termed  the  inverse-velocity  gradient  [20],  In  the  apex  region,  each  jet 
element  has  a  higher  velocity  than  the  one  ahead  of  it,  causing  a  "piling-up" 
of  jet  mass.  This  piled  up  mass  forms  the  Jet  tip  as  verified  experimentally 
[21].  For  conventional  conical  charges,  the  first  30  to  *10$  of  the  liner  from 
the  theoretical  apex  forms  the  tip  of  the  Jet.  The  tip  has  a  much  larger 
radius  than  the  rest  of  the  jet. 


Figure  9  shows  the  inverse-velocity  gradient  with  respect  to  X,  the  liner 
element  position.  The  velocity  of  the  jet  tip  is  given  as 


V. 

jo 


V, 


TIP 


djn  •  j  v 
3  dX 

dX 


dx 


> 


based  on  the  conservation  of  momentum. 


(33) 


Recently,  Hirsch  [22]  extended  the  PER  model  to  include  non-steady 

effects  at  the  collision  point  (or  formation  point)  by  assuming  that  relative 

to  this  point,  the  jet  and  slug  velocities,  V.  and  V_,  are  not  equal  to  the 

J  3 

incoming  liner  velocity  V  as  the  PER  theory  assumed  from  the  Bernoulli 
equation. 

Other  studies  involve  modifications  or  extentions  of  the  PER  theory,  such 
as  Perez,  et  al  [233  who  modified  PER  to  account  for  two-dimensional  effects 
of  the  flow  near  the  collapse  axis,  and  Hirsch  [243  .  Leidel  [253  provides  an 
overview  of  the  PER  model  and  analyzes  annular  cutting  charges,  or  cookie 
cutter  warheads  as  first  conceived  by  Glass,  et  al  [263.  Harrison  [27,  283 
compares  the  simple  one-dimensional  models  to  hydrocodes  and  experimental 
data.  These  methods  are  also  related  to  warhead  design. 

Other  shaped  charge  publications  of  interest  include  Hornemann  and  Senf 
[293  on  wide  angle  cones.  Van  Thiel  and  Levatin  [303,  Trinks  [31]  for  his 
shaped  charge  discussion  and  historical  inferences,  Rostoker  [323,  Simon  and 
DiPersio  [333,  and  the  USSR  studies  by  Novikov  [34]  and  Titov  [353.  These 
papers  relate  to  jet  formation  and  collapse  and  related  topics  pursuing  the 
concepts  we  have  reviewed  so  far.  An  excellent  Soviet  paper  was  published  by 
Ivanova  and  Rozantseva  [363  which  is  very  long  and  very  informative.  It 
provides  a  detailed  discussion  and  derivation  of  the  PER  theory  and  simple 
penetration  models  including  experimental  data  and  numerous  flash  radiographs 
or  x-ray  photographs  of  shaped  charge  Jets.  The  USSR  models  are  discussed  in 
more  detail  later  in  the  text. 


IX-20 


Excellent  flash  radiography  data  and  analysis  thereof  is  given  by 
Breidenbach  [37]  and  by  Heine-Geldern  and  Pugh  [38].  Classical  papers  on  jets 
and  jet  formation  are  given  by  Harlow  and  Pracht  [39]  and  on  collapsing 
cylinders  by  Koski,  et  al  [40]. 

Zernow  and  Simon  [41].  investigated  the  behavior  of  jets  from  several 
different  shaped  charge  materials.  Aseltine  [42]  studied  the  effect  of 
several  different  warhead  asymmetries  on  the  behavior  and  performance  of 
shaped  charge  jets.  Merendlno,  et  al  [43]  altered  a  shaped  charge  liner, 
i.e..  inhibited  the  collapse  of  the  liner  to  obtain  a  massive  hypervelocity 
pellet  in  lieu  of  a  jet.  Other  studies  regarding  shaped-charge  analyses  are 
given  in  the  recent  papers  by  Robinson  [44.  45]  and  Lee  [46]. 

Kolsky.  et  al  [47]  studied  the  formation  and  penetration  of  conical 
shaped  charge  liners.  Kolsky  [483  continued  his  studies  regarding  the 
formation  and  penetration  of  hemispherical  shaped  charge  liners.  The 
formation  and  collapse  of  hemispherical  liners  were  also  studied  via 
hydrocodes  by  Kiwan  and  Arbuckle  [49],  Arbuckle,  et  al  [50],  Chou,  et  al  [51] 
Aseltine,  et  al  [52],  Chou,  et  al  [533.  and  Lacetera  and  Walters  [54], 

Studies  on  the  formation  of  hemispherical  shaped  charge  liners  were 
presented  by  Lee  [46]  Singh  [55],  Shepherd  [56]  and  Grace  [57].  A  closed  form 
analytical  expression  (analogous  to  the  PER  theory  for  conical  liners)  for  the 
collapse  and  formation  of  hemispherical  shaped-charge  liners  does  not 
currently  exist. 
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X.  THE  VISCO-PLASTIC  THEORY 


Other  studies  involving  the  collapse  and  formation  of  shaped  charge 
liners  were  advanced  by  Soviet  researchers.  Returning  to  the  collapse  of 
conical  liners  or  wedges,  the  USSR  used  the  results  of  Birkhof f ,  or  perhaps 
Lavrent'ev,  as  the  basis  for  their  studies,  Lavrent'ev  [1]  having  claimed 
prior  discovery.  In  any  case,  the  results  of  Birkhof f  or  Lavrent'ev  were 
extended  to  include  visco-plastic  effects  by  Godunov,  et  al  [2]. 

The  visco-plastic  model  follows  from  References  [1  -  7]  and  a  detailed 
derivation  and  discussion  of  this  model  is  given  by  Walters  [8].  The  basic 
viscoplastic  theory  [1  -  7]  does  not  include  transient  effects.  Typically,  in 
the  axisymmetric  hydrocode  models  used  at  the  BRL  and  other  agencies, 
compressible  flow  is  assumed  and  the  material  constitutive  relationships  are 
based  on  elastic-perfectly  plastic,  work  hardening  models  or  more 
sophisticated  models,  e.g.  [9,  10].  The  USSR  analytical  models  typically 
assume  an  Incompressible  flow,  but  use  a  rate-dependent,  visco-plastic 
material  constitutive  equation  [2,  3,  11]*  This  equation  is  usually  of  the 

form  o=o  +  vt,  where  the  stress  o  is  related  to  the  strain  rate  £,  by  the 

y 

yield  stress  °y  and  the  constant  dynamic  viscosity  coefficient,  b. 

Thus,  the  visco-plastic  models  require  a  knowledge  of  the  dynamic 
viscosity  coefficient  and  many  USSR  investigators  have  deduced  viscosity 
coefficients  from  experimental  measurements  under  shock  loading  conditions. 
Walters  [11]  summarizes  many  of  the  experimental  viscosity  values.  Additional 
data  related  to  viscosity  measurements  and  the  collapse  of  metallic  bodies  by 
high  explosives  are  given  by  References  [12  -  34].  Material  viscosity  values 
have  also  been  experimentally  determined  in  the  US  [11,  35],  under  shock 
loading  conditions.  The  viscosity  values  deduced  from  various  experimental 
measurements  depend  on  many  parameters,  primarily  strain  rate,  pressure  and 
temperature.  Thus,  the  viscosity  coefficient  is  not  constant.  In  fact,  for 

5 

solid  metals,  the  dynamic  viscosity  may  range  from  typically  10  to  10  Pa-sec 

2  6 

(10  -  10  poise)  depending  on  the  strain  rate,  pressure  and  temperature  and 

the  experimental  method  used  to  measure  b.  Usually,  the  US  measured  viscosity 
values  under  shock  loading  conditions  may  be  as  much  as  two 
orders-of-magnitude  lower  than  the  USSR  measured  viscosity  values  [11  -  35]. 
The  uncertainty  in  the  US  measured  b  values  is  +_  a  factor  of  ten.  The  USSR 
methods  may  even  have  a  greater  amount  of  uncertainty.  Experimental 
determination  of  b  in  the  US,  as  performed  by  Chhabildas  and  Asay  [35], 
involves  measurement  of  the  shock  rise  time  via  velocity  interferometry.  The 

maximum  stress  is  taken  to  be  proportional  to  the  Btrain  rate,  o'vpt,  and  the 
mflYimum  stress  is  determined  from  the  Rayleigh  line  and  the  shock  Hugoniot 
pressures.  The  strain  rate  is  calculated  as  the  strain  induced  by  the  shock 
divided  by  the  shock  rise  time,  which  is  resolution  limited.  The  time 
resolution  is  limited  to  1  to  3  nanoseconds  [35]. 

The  USSR  shaped-charge  liner  collapse  model  includes  the  influence  of 
viscosity  on  the  jet  formation  process.  A  Jet  formation  criterion  is  derived 
which  is  based  not  on  taking  account  of  the  compressibility  effects  or  shock 


X-l 


effects  as  in  References  [36]  or  [37],  or  on  consideration  of  a  critical  Mach 
number  as  in  Reference  [38],  but  on  taking  into  account  the  viscous  properties 
of  metals. 

Following  Godunov,  et  al  [2],  an  appropriate  method  is  used  to  estimate 
the  effect  of  viscosity  in  the  plane  problem  (two-dimensional)  of  Jet 
collisions.  The  fluid  is  incompressible,  the  motion  irrotational  and  steady 
state,  and  the  coefficient  of  viscosity  is  constant.  Then  the  solutions  of 
Euler  equations  automatically  satisfy  the  Navier-Stokes  equations  of  motion 
and  the  difference  between  the  problems  of  ideal  and  viscous  Jet  collisions  is 
the  conditions  on  the  free  surface  of  the  Jet.  The  flow  fields  are  taken  to 
be  Newtonian.  For  the  flow  fields  to  agree  in  the  ideal  and  viscous  flow 
problems,  some  forces  must  be  applied  to  the  free  surface  in  the  viscous  flow 
case  [2],  Consideration  of  the  influence  of  these  forces  on  the  flow 
resulting  from  the  collision  between  Jets  of  an  ideal  fluid  will  provide  a 
measure  of  the  influence  of  viscosity  in  the  Jet  formation  problem  [2]. 

The  formation  of  a  reverse  Jet  (slug)  is  possible  only  if  the  horizontal 
component  of  the  viscous  force,  acting  on  the  free  surfaces  of  the  reverse 
jet,  is  less  that  the  force  resulting  from  the  symmetric  collision  of  two 
plane  fluid  Jets  [2].  This  inequality  allows  one  to  define  a  Reynolds  number 
and  a  critical  Reynolds  number  for  jet  formation. 

The  visco-plastic  shaped-charge  jet  oollapse  model  derived  by  Godunov  [2] 
relaxes  to  the  classical  US  models  reviewed  earlier  when  the  dynamic  viscosity 
is  zero.  Basically,  the  visco-plastic  model  yields  a  lower  Jet  velocity  and  a 
lower  flow  velocity  than  the  US  models.  However,  the  visco-plastic  model 
predicts  a  higher  slug  velocity  than  the  US  models.  The  USSR  (visco-plastic) 
criterion  to  form  a  coherent  jet  is  that  the  Reynolds  number, 


Re  = 


tusin  6 


v  (1  -  sing) 


(1) 


be  greater  than  two.  The  wall  thickness  is  denoted  by  t,  v  is  the  kinematic 
viscosity,  2 6  is  the  collision  angle,  and  u  is  the  inviscid  flow  velocity  as 
given  by  Defourneaux  or  the  PER  model  for  example.  The  Reynolds  number 
greater  than  two  criterion  is  based  on  experimental  observations  [2,  11]. 
Detailed  USSR  formulae  for  the  Jet  strain  rate,  jet  velocity,  viscid  flow 
velocity,  and  slug  velocity,  are  given  in  References  [2,  11].  Also,  by 
setting  Re  =  2  in  Equation  1 ,  an  expression  for  the  critical  flow  velocity  as 
a  function  of  (5,  t  and  v  can  be  obtained.  If  the  critical  flow  velocity  is 
greater  than  the  flow  velocity,  a  coherent  jet  will  not  form.  Simonov  [6]  and 
others  [11  -  14,  18,  20,  21,  24,  27,  28,  31,  32]  define  collapse  angles  for 
the  transition  from  the  jet  regime  to  the  jetless  regime  for  explosive  bonding 
applications.  These  USSR  investigators  defined  collapse  angles  below  which 
only  incoherent  jets  will  form  and  collapse  angles  below  which  no  jet  at  all 
will  form.  Walsh,  et  al  [36]  served  as  the  basis  for  many  of  the  USSR 
studies. 
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For  the  transient  flow  (PER)  case,  the  stagnation  point  velocity  (V  )  is 
given  a3  ^ 


V 


c 


sin(g  -  q)-sin(g  -  a  -  <ft) 
sing 


(2) 


where  D  is  the  speed  of  the  detonation  wave,  2a  is  the  conical  liner  apex 
a 

angle  and  $  =  26  is  the  plate  bending  angle.  Equation  2  is  used  in  both  the 
PER  and  visco-plastic  models.  The  flow  velocity  is  given  by 


si 


sin(a  +  jQ-sina 


sing 


i 


(3) 


in  the  PER  model.  The  visco-plastic  flow  velocity  is 


U2  =  u  Cl  -  2/Re) ^ 


(*0 


from  Reference  [2]  for  Re  >  2.  Also,  U2  =  U  when  the  dynamic  viscosity,  y  ,  is 
zero. 


In  the  PER  model,  the  jet  velocity  and  slug  velocity  are,  respectively, 

V,  =  V  +  U,  (5) 

j  c 

and 

V  =  V  -  U.  (6) 

s  c 

In  the  visco-plastic  model,  the  jet  and  slug  velocities  are 

Vj  *  ’o  *  V  (7) 

and 


V  =  V  -  U_. 
s  c  2 

In  both  models, 


(8) 


V,  +  V  =  2V  . 

J  s  c 

For  steady  state  collapse,  analogous  formulae  are  available  from  Godunov 
and  Birkhoff. 


For  both  the  steady  state  and  transient  version  of  the  visco-plastic 
theory,  the  US  theory  can  be  obtained  by  setting  the  dynamic  viscosity,  y,  or 
the  kinematic  viscosity,  v  =  u/p  equal  to  zero.  Walters  and  Harrison  [8] 
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discuss  the  PER  and  the  visco-plastic  models  and  comparison  is  made  to  the 
experimental  data  obtained  by  Allison  and  Vitali  [39]. 

The  strain  rates  calculated  by  the  visco-plastic  model  [2]  are  of  the 

order  of  10  ^/s  and  show  little  variation  with  liner  element  position.  Chou 
and  Caleone  [MO]  calculated  strain  rates  for  a  81  mm  diameter  precision  charge 

with  a  copper  conical  liner  with  a  M2°  apex  angle  and  a  1 .9  mm  wall  thickness 
and  found  a  larger  variation  of  strain  rate  with  a  liner  element  position 
including  liner  elements  in  both  compression  and  tension.  The  peak  strain 

5 

rate  was  10  /s.  Experimental  measurements  of  the  strain  rate  of  collapsing 
liners  are  not  available,  but  Bauer  and  Bless  [Ml]  measured  the  strain  rate  of 

ti 

exploding  copper  tubes  to  be  about  10  /s  and  Walters  deduced  the  strain  rates 

M 

of  stretching  Jets  to  be  of  the  order  of  10  /s. 

Comments  on  .Jet,  Coherency. 

The  visco-plastic  coherent  Jet  or  incoherent  jet  criterion  is  based  on  a 
critical  Reynolds  number  being  greater  than  2  [2],  or  some  other  number  [MM]. 
The  US  criterion  is  inviscid,  but  compressibility  effects  are  taken  into 
consideration  in  the  determination  of  the  cohesiveness  of  the  jet. 

Walsh,  et  al  [36]  studied  the  compressibility  effect  in  the  jet  formation 

process  as  related  to  the  oblique  impact  of  explosively  driven  plates.  The 

symmetric  plate  collision  as  viewed  from  a  moving  coordinate  system,  reduces 

to  a  flow  field  analogous  to  two  impinging  streams.  Walsh,  et  al  [36] 

concluded  that  Jetting  always  occurs  if  the  fluid  is  incompressible  or  if  the 

collision  velocity  in  the  moving  coordinate  system  is  subsonic.  For 

supersonic  flow,  Jetting  always  occurs  if  P>  where  0  is  the  critical 

c  c 

turning  angle  for  an  attached  oblique  shock  wave  at  the  collision  velocity. 
Walsh,  et  al  [36]  did  not  address  jet  cohesiveness  per  se,  since  their  primary 
interest  was  explosive  bonding.  Cowan  and  Holtzmann  [M5 ]  present  additional 
flow  criteria  for  explosive  bonding. 

In  the  studies  of  Chou,  et  al  [37,  M6],  Jetting  criteria  for  plane 
axisymmetric  cases  are  presented  along  with  a  measure  of  the  jet  quality. 

From  References  [37]  and  [M6],  Chou  states: 

1.  For  subsonic  collisions  (or  the  collision  velocity  V  <  C,  the 
material  bulk  speed  of  sound)  a  solid  coherent  jet  always  forms; 

2.  For  supersonic  collisions  (V  >  C)  Jetting  occurs  if  P  >  /3  ,  but  the 

c 

jet  is  not  coherent.  The  angle  PQ  is  the  maximum  angle  that  an  attached  shock 
wave  can  form  at  a  prescribed  supersonic  velocity,  V;  and 

3.  For  supersonic  collisions  (V  >  C)  but  P  <  P ,  a  Jet  will  not  be 

formed.  0 

For  shaped  charge  applications,  the  major  criterion  for  Jetting  is  that 
the  formation  process  be  subsonic.  Otherwise,  the  jet  will  be  incoherent  and 
spread  out  radially. 


The  speed  of  sound  referred  to  is  usually  taken  to  be  the  bulk  speed  of 
sound  of  the  material  as  opposed  to  a  longitudinal  or  transverse  shear  speed 
of  sound.  The  appropriate  speed  of  sound ,  as  well  as  its  exact  value  under 
the  extreme  pressures  and  temperatures  encountered  in  the  formation  region  is 
not  well  known.  Also,  the  actual  flow  field  is  compressible  and 
two-dimensional  (at  least)  and  hence  the  subsonic  collision  condition  is 
useful  only  as  a  general  principle.  In  practice,  it  has  been  observed  that 
the  value  of  the  critical  Mach  number  based  on  the  static  bulk  speed  of  sound 
is  about  1.2  for  conical  copper  liners.  Harrison  [38.  *17].  and  Walters  and 
Harrison  [8].  In  other  words,  the  liner  velocity  divided  by  the  bulk  speed  of 
sound,  i.e.,  the  Mach  number,  could  be  as  high  as  1.2  and  still  achieve  a 
cohesive  jet  for  a  conical,  copper  liner.  A  flow  velocity  of  4.8  km/s  (1.23 
times  the  bulk  speed  of  sound  of  copper)  is  the  calculated  flow  velocity  above 

which  20  copper  liners  have  been  observed  to  have  incoherent  Jet  tips  [38]. 
Figure  1  shows  a  flash  radiography  photograph  of  a  copper,  coherent  and 
incoherent  jet. 
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Figure  1.  Radiographs  of  Jets  From  Two  Typical  Conical  Charges 

(a)  40°  Copper-lined  Charge,  Subsonic  Collision, 
Coherent  Jet,  and 

(b)  20°  Copper-lined  Charge,  Supersonic  Collision, 
Incoherent  Jet. 
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XI.  EXPLOSIVE  WELDING,  SHAPING  AND  FORMING 

As  mentioned  in  the  first  section,  we  will  briefly  address  the  concepts 
of  explosive  welding,  explosive  bonding  and  explosive  forming. 

The  explosive  working  of  metals  and  concepts  related  to  explosive 
effects,  standoff  systems,  contact  systems,  shock  physics  and  material  effects 
are  presented  in  an  excellent  introductory  book  by  Rinehart  and  Pearson  til. 

An  explosive  forming  procedure  is  shown  in  Figure  1 .  A  die  assembly 
holding  the  work  piece  is  positioned  in  a  tank  of  water.  A  vacumm  is  drawn 
below  the  work  piece  and  an  explosive  charge  is  detonated  in  the  water.  The 

work  piece  then  assumes  the  shape  of  the  die  (if  all  goes  well). 

A  method  for  bulging  cylindrical  tubing  is  shown  in  Figure  2.  The 
explosive  charge  is  sealed  in  plastic  and  positioned  in  the  cylindrical  work 
piece  filled  with  water.  The  ends  are  sealed  to  add  confinement  and  reduce 
water  splash.  The  charge  is  detonated  within  the  cylinder  held  by  a  split  die 
assembly.  Figure  2  shows  the  bulging  cylinder  produced. 

In  explosive  welding  or  bonding  applications,  two  or  more  plates  are 
driven  together  at  high  velocities  and  surface  jetting  occurs.  This  results 
in  a  bond  between  the  two  driven  metals  and  can  result  in  a  wave  interface  of 
high  strength  between  the  welded  materials.  If  the  plate  impact  velocities 
are  too  high  (relative  to  the  material  properties),  the  interface  material 
will  jet  as  a  shaped-charge  and  the  bonding  will  be  poor.  If  the  impact 
velocities  are  too  low,  no  interaction  will  occur  at  the  interface.  The 
velocity  must  be  such  that  plastic  deformation  occurs  at  the  interface  and 
local  surface  Jetting  occurs.  The  requirements  for  welding  to  occur  are  [2): 

a.  the  existence  of  a  jet  at  the  interface,  and 

b.  an  increase  in  pressure,  associated  with  the  rapid  dissipation  of 
kinetic  energy,  to  a  sufficient  level  for  a  sufficient  time  to  achieve  stable 
interatomic  bonds.  The  pressure  is  determined  by  the  impact  velocity  and  the 
time  available  for  bonding  is  determined  by  the  velocity  of  the  collision 
point . 


The  final  geometry  of  the  explosive-metal  system  will  vary  depending  on 
the  geometry  and  materials  of  the  parts  to  be  welded.  The  nature  of  the 
interface  between  the  welded  parts  will  vary  depending  on  the  materials  welded 
and  the  nature  of  the  explosive-metal  interaction  system. 

Figure  3  shows  a  surface  cladding  test  arrangement.  The  top  plate,  or 
flyer  plate,  only  is  driven  by  the  explosive  force  in  order  to  weld  a  thin 
sheet  of  material  onto  a  heavy  sheet.  This  technique  is  useful  in  plating  a 
good  structural  material  with  a  thin  cladding  to  protect  it  from  a  corrosive 
or  hazardous  environment. 

Explosive  welding  is  similiar  to  explosive  cladding,  but  two  flyer  plates 
are  used  to  form  an  explosively  welded  final  product.  Figure  4  depicts  an 
explosively  welded  sandwich  or  a  three-layer  weld. 

Other  examples  of  the  explosive  working  of  metals  as  well  as  the  details 
associated  with  these  processes  can  be  found  in  Reference  [1], 
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Blazynski  [2]  provides  a  recent  collection  of  papers  dealing  with  many 
aspects  of  the  explosive  working  of  metals.  Pearson,  in  Chapter  1  of 
Reference  [2],  describes  the  explosive  hardening  of  austenitic  manganese 
steels  subjected  to  severe  impact  and  abrasion.  These  steels  are  used  in 
railroad  frogs,  rock  crusher  jaws,  grinding  mills  and  similar  devices.  To 
explosively  harden  a  steel,  a  thin  layer  of  explosive  is  detonated  in  contact 
with  the  surface  to  be  hardened,  usually  at  grazing  incidence.  The 
propagation  of  the  shock  through  the  material  increases  the  metal  hardness 
(BHN)  throughout  the  material.  This  also  increases  the  yield  and  tensile 
strength  of  the  metal,  making  it  more  resistant  to  impact  wear  and  with  less 
tendency  to  deform.  Duvall  [3]  discusses  shock  wave  propagation  in  metals. 

Pearson,  in  Reference  [2],  discusses  explosive  compaction  where  explosive 
loads  are  applied  either  directly  or  through  a  loading  system  to  compact 
powders.  Two  methods  involve  the  use  of  gun  powder  cartridges  and  the  use  of 
explosives  to  drive  pistons  into  the  powder  compaction  chamber.  Explosive 
compaction  can  produce  parts  from  powders  that  cannot  be  produced  by 
conventional  pressing.  Improved  material  properties  with  high  press  densities 
can  be  produced.  Explosive  compaction  is  extremely  valuable  in  the  powder 
metallurgy  field. 

Pearson,  Chapter  1  and  El-Sobky,  Chapter  6  of  Reference  [2],  discuss 
explosive  welding.  We  have  briefly  touched  on  this  area  and  we  will  conclude 
our  discussion  with  some  of  the  major  points  of  Reference  [23. 

Explosive  welding  is  a  solid  phase  welding  process  in  which  high 
explosives  are  used  to  Join  the  weld  surfaces  in  a  high  velocity  collision. 
This  collision  produces  severe,  localized  plastic  flow  at  the  interface 
between  the  two  surfaces.  Explosive  welding  is  used  to  weld  metal 
combinations  many  of  which  cannot  be  welded  by  conventional  means. 

The  weld  is  usually  formed  by  oblique  impact  of  the  plate  surfaces  with 
the  weld,  progressing  from  the  apex  along  the  collision  interface  by  a  process 
similar  to  the  collapse  of  a  shaped-charge  liner.  The  parameters  involved  in 
explosive  welding  include  the  physical  and  mechanical  properties  of  the  metals 
to  be  welded,  the  type  and  amount  of  explosive  used,  the  mode  of  initiation, 
the  initial  geometry  of  the  weld  operation  and  the  type  and  geometry  of  buffer 
sheets  (if  used)  between  the  metal  and  the  explosive.  These  parameters 
influence  the  collision  angle,  the  impact  velocity  and  the  collision  point 
velocity. 

It  is  generally  agreed  that  jet  formation  at  the  collision  point  is 
essential  for  welding,  e.g.,  El-Sobky,  Chapter  6  of  Reference  [23.  This  is 
because,  at  least  in  part.  Jetting  produces  chemically  clean  surfaces  which 
allow  interatomic  bonding  under  the  high  pressure  conditions. 

The  jet-no- Jet  criterion,  or  the  limiting  condition  for  jet  formation, 
discussed  earlier,  is  based  on  the  studies  of  Cowan  and  Holtzman  [43,  Walsh, 
et  al  [53,  Chou,  et  al  [63  and  Harlow  and  Pracht  [73.  These  studies  relate  to 
shock  formation  and  critical  collapse  angles  and  are  also  applied  to 
shaped-charge  jets.  The  mathematical  description  for  the  explosive  welding  of 
oblique  plates  involves  relationships  similar  to  those  used  in  jet  formation 
theory  and  in  the  Gurney  type  equations.  In  addition  to  Gurney,  impact 
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velocity  equations  are  given  by  Deribas,  et  al  [8]  and  others  as  reviewed  by 
El-Sobky ,  Chapter  6  of  Reference  [2].  El-Sobky  reviews  the  jet  formation 
criterion  as  well  as  minimum  impact  pressures  and  standoff  distance 
requirements  which  must  be  met.  El-Sobky  notes  that  the  collision  velocity 
and  plate  impact  velocity  should  be  less  than  the  bulk  speed  of  sound  of 
either  welding  component.  However,  Wylie,  et  al  [9]  suggests  that  the  bulk 
speed  of  sound  may  be  exceeded  by  as  much  as  25?  with  satisfactory  welds.  It 
is  interesting  that  this  coincides  with  the  jet-no- jet  criterion  observed  by 
Harrison,  et  al  [10,11,12]  for  copper  jets,  i.e.,  Mach  number  1.23.  Also,  the 
visco-plastic  critical  Reynolds  number  criterion,  discussed  earlier,  was 
derived  from  explosive  welding  experiments  [13,14]. 

A  major  area  of  study  in  explosive  welding  is  the  interfacial  periodic 
deformation  or  the  interfacial  waves.  Often  a  regular  wavy  interface  pattern 
appears,  but  the  nature  of  the  interfacial  waves  vary  depending  on  the 
materials  and  system  geometry  used. 

When  regular  waves  are  produced,  they  are  very  similar  to  the  surface 
deformations  found  on  solid  surfaces  subjected  to  erosive  action  of  high  speed 
liquid  films,  jet  drags,  e.g.,  in  pump  and  steam  turbine  blades,  and  at  the 
air-sea  interface,  e.g.,  ocean  waves.  A  review  of  the  mechanisms  for  wave 
formation  are  given  by  El-Sobky  in  Chapter  6  of  Reference  [2], 

A  good  general  discussion  of  explosive  welding  is  given  in  Reference 
[15].  Jetting  phenomena,  the  weld  interface,  and  the  explosive  load  and 
standoff  distance  required  for  various  weld  combinations  are  discussed.  Other 
excellent  sources  of  information  include  the  USSR  sources,  given  in  the 
discussion  on  the  visco-plastic  Jet  formation  model  and  viscosity  measurements 
[13»14].  Schroeder  [ 1 6 ]  recently  edited  a  text  on  explosive  welding  and 
related  areas.  Deribas  [17]  authored  a  classic  text  on  explosive  welding  and 
work-hardening.  Duvall  [3]  presents  an  excellent  supplementary  discussion  on 
shock  propagation  and  transmission  through  solids.  Sewell  [18]  also  discusses 
explosive  welding  regarding  jet-no- jet  criterion. 

This  concludes  the  discussion  of  explosive  working  of  metals.  The 
information  provided  in  this  section  is  not  complete,  nor  was  it  intended  to 
be.  The  intent  was  merely  to  provide  an  overview  of  explosive  welding  and 
related  fields.  Our  goal  was  to  provide  a  brief  introduction  and  a  few 
fundamental  references  into  this  field  and  to  note  the  similarity  to  shaped- 
charge  Jetting  phenomenon. 
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Figure  1.  Explosive  Forming, 
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Figure  2.  Bulging  Cylinder  Tube. 
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